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Evaluations  of  Duck  Habitat  and  Estimation  of  Duck  Population 
Sizes  with  a  Remote-Sensing-Based  System 
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National  Biological  Ser\nce 
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Abstract.  During  1987-90,  we  used  high-altitude  photography,  aerial  videography,  counts,  and 
models  to  estimate  sizes  of  breeding  populations  of  dabbling  ducks  (Anatinae)  and  duck  production 
and  to  identify  duck  habitat  on  U.S.  Fish  and  Wildlife  Service  land  and  easements  and  on  private  land 
in  the  prairie  pothole  region  of  the  United  States.  The  study  area  contained  about  3.1  million  wetland 
basins  (28,490  km^).  Wetland  area  (ha  per  km^)  was  highest  on  service-owned  land;  wetland-basin 
density  was  greatest  on  service  easements.  Temporary  and  seasonal  wetlands  were  underrepresented 
and  lakes  were  overrepresented  on  service-owned  land.  Seventy-eight  percent  of  all  basins  were  less 
than  0.41  ha.  Cropland  dominated  private  land.  Pond  density  decreased  from  4.4/km^  in  1987  to 
3.4/km^  in  1990  and  pond  area,  from  7,2  ha/km^  to  2.7  ha/km^.  The  density  of  the  blue-winged  teal 
was  greatest  (3.4  pairs/km^)  and  was  followed  in  magnitude  by  those  of  the  mallard  (2.1  pairs/km  ), 
the  gadwall  (1 .8  pairs/km^),  the  northern  pintail  (0.8  pairs/km^),  and  the  redhead  (0.8  pairs/km^).  Duck 
density  was  consistently  highest  on  service-owned  land.  The  decline  of  breeding-population  sizes  in 
1987-90  closely  corresponded  to  losses  of  pond  numbers  and  pond  area.  The  density  of  breeding  pairs 
per  pond  was  inversely  related  to  pond  density,  suggesting  that  breeding  ducks  tended  to  concentrate 
on  the  remaining  ponds  as  drought  intensified.  The  production  of  recruits  followed  the  same  pattern 
as  breeding-population  sizes.  We  estimated  that  2.5%  of  the  ducklings  hatched  on  service-owned  land, 
which  was  1.3%  of  the  study  area;  19.6%  hatched  on  service  easements,  which  were  14.2%  of  the 
study  area;  and  77.9%  hatched  on  private  land,  which  was  84.6%  of  the  study  area.  Various  sources 
of  bias  and  sampling  error  and  improvements  to  the  system  are  discussed. 

Key  words:  Ducks,  wetlands,  population,  recruitment,  management,  models,  videography,  drought. 


The  U.S.  Fish  and  Wildlife  Service  has  long  been  in¬ 
volved  in  monitoring  wetland  conditions  and  duck  popula¬ 
tions.  Cooperative  breeding-ground  surveys  by  the  United 
States  and  Canada  (Martin  et  al.  1979;  Reynolds  1 987)  have 
provided  data  essential  to  management  of  continental  water- 
fowl  populations.  However,  these  surveys  were  designed  for 
obtaining  data  that  are  used  primarily  for  setting  annual 
hunting  regulations  and  were  not  intended  for  obtaining  data 
for  the  management  of  national  wildlife  refuges  and  wetland 
management  districts  or  for  evaluating  differences  of  habitat 
use  by  ducks  among  landownership  classes. 

To  meet  the  unfilled  need  for  data  from  local  areas, 
Hammond  (1969)  developed  a  system  for  measuring  duck 
use  and  production  on  refuges  and  in  waterfowl  production 
areas  in  the  prairie  pothole  region  of  the  United  States.  That 


system  was  frequently  modified  to  meet  the  needs  and 
resources  of  individual  managers,  and  some  of  the  parame¬ 
ters  that  were  used  in  the  procedure  were  later  shown  to  be 
erroneous.  In  addition,  the  Hammond  system  did  not  allow 
comparisons  of  duck  use  and  production  among  lands 
owned  by  the  U.S.  Fish  and  Wildlife  Service,  easements  by 
the  service,  and  private  lands.  This  is  a  critical  need  by 
managers  in  the  prairie  pothole  region.  With  our  system,  we 
attempted  to  overcome  some  of  the  shortcomings  of  the 
previous  system. 

Here  we  describe  a  remote-sensing-based  system  for 
estimating  the  number  and  area  of  ponds,  the  sizes  of  breed¬ 
ing  duck  populations,  and  the  number  of  young  ducks  re¬ 
cruited  to  the  population  each  fall  in  the  prairie  pothole 
region  of  Minnesota,  North  Dakota,  South  Dakota,  and 
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Montana.  The  area  of  upland  habitat,  number  of  ponds,  and 
areas  of  ponds  in  18  waterfowl  management  districts  (all 
districts  in  the  Dakotas  and  in  Minnesota  but  only  the  most 
eastern  districts  in  Montana)  were  estimated. 

Our  goals  were  to  design  a  system  that  provides  consis¬ 
tent  estimates  among  areas  and  years,  to  maximize  the  use 
of  existing  data,  to  rapidly  execute  complex  procedures 
with  microcomputers,  and  to  develop  a  database  that  docu¬ 
ments  changes  in  the  estimated  parameters. 

Methods 

Definitions 

Definitions  of  the  terms  wetland,  wetland  area,  wetland 
basin,  and  pond  follow  Cowardin  (1982).  We  define  four 
essential  terms  for  describing  nesting  ducks:  A  successful 
clutch  is  a  clutch  of  eggs  of  which  one  or  more  eggs  hatch 
even  if  one  or  more  dead  young  are  found  at  the  nest  site 
(see  successful,  Klett  et  al.  1986).  A  successful  nesting 
attempt  is  a  hen’s  attempt  to  nest  when  that  attempt  results 
in  a  successful  clutch.  Clutch  success  is  the  probability  that 
a  clutch  will  be  successful  and  is  directly  equivalent  to  nest 
success  as  defined  by  Klett  et  al.  (1986).  Hen  success  is 
the  probability  that  a  hen  will  make  at  least  one  successful 
nesting  attempt  during  a  single  breeding  season. 

Design 

The  remote-sensing-based  system  differs  from  previous 
systems  in  concept,  scope,  and  objectives.  For  instance, 
there  is  a  relation  between  the  amount  of  wetland  habitat  in 
an  area  and  the  number  of  ducks  expected  to  settle  there  in 
spring.  Johnson  and  Grier  (1988)  presented  a  rigorous  dis¬ 
cussion  of  the  relation  between  numbers  of  ponds  and 
numbers  of  ducks  observed  during  the  cooperative  breed¬ 
ing-ground  surveys.  A  less  obvious  relation  is  between  the 
amount  and  quality  of  nesting  habitat  and  the  productivity 
of  ducks  (Klett  et  al.  1 988).  Direct  measures  of  the  breeding 
population  and  production  are  difficult  and  expensive. 
Measurements  of  the  amount  and  type  of  habitat  by  remote 
sensing  are  relatively  simple  and  inexpensive.  Therefore, 
models  that  relate  duck  numbers  and  production  to  habitat 
should  increase  the  precision  of  surveys  without  greatly 
increasing  the  cost. 

Johnson  et  al.  (1987)  developed  a  mallard  productivity 
model.  That  model  and  habitat  data  from  a  large  sample 
{n  =  422)  of  10.4-km^  plots  and  nest-survival  estimates  in 
those  habitats  (Klett  et  al.  1988)  were  used  by  Cowardin 
et  al.  (1988)  to  simulate  results  from  various  management 
scenarios.  These  plots  and  the  data  files  from  them  were 
used  as  the  basis  for  the  remote-sensing-based  system. 
From  1982  to  1986,  preliminary  compilation  of  data 
and  tests  of  proposed  techniques  were  conducted  in  the 


Arrowwood  Waterfowl  Management  District.  The  study 
included  breeding-pair  counts  on  10.4-km^  plots  and 
building  of  baseline  regression  equations  for  estimating 
duck  numbers  from  pond  data.  These  regressions  were 
specific  to  areas  and  years  when  data  were  available.  They 
were  later  modified  to  account  for  annual  and  regional 
variations.  We  also  assessed  the  adequacy  of  the  regres¬ 
sion  equations  for  estimating  duck  numbers  and  evaluated 
video  cameras. 

Sample  Universe 

The  remote-sensing-based  system  was  applied  in  the 
prairie  pothole  region  of  the  United  States  in  Minnesota, 
North  Dakota,  South  Dakota,  and  Montana.  This  area  of 
glacial  landscape  is  bounded  on  the  east  by  forest  land  in 
Minnesota,  on  the  south  and  west  by  the  limit  of  glaciation 
in  the  Dakotas  and  Montana,  and  on  the  north  by  the 
Canadian  border  (Fig.  1).  We  approximated  the  bounda¬ 
ries  by  transferring  boundaries  presented  by  Hammond 
(1965)  and  Mann  (1974)  to  1:500,000  U.S.  Geological 
Survey  (USGS)  maps  with  the  constraint  that  boundaries 
must  follow  townships,  which  were  used  as  a  basis  for 
stratifying  sampling  units. 

Sample  Design 

The  sampling  units  for  habitat  data  were  10.4-km^ 
plots.  The  plot  size  was  chosen  to  approximate  the  home- 
range  size  of  a  breeding  mallard  pair  (Cowardin  et  al. 
1988).  By  1990,  the  sample  of  335  plots  in  1987  was 
increased  to  443  plots  (Table  1). 

The  original  sample  was  a  stratified  random  sample  of 
500  plots  drawn  from  a  universe  that  represented  the  United 
States  portion  of  the  prairie  pothole  region  (Cowardin  et  al. 
1988;  Fig.  2).  The  universe  was  divided  into  93.2  km^ 
townships  based  on  1:500,000  state  maps.  We  defined  3 
landownership  classes:  land  owned  in  fee  by  the  U.S.  Fish 
and  Wildlife  Service  {service-owned  land);  easements  con¬ 
sisted  of  tracts  for  which  the  service  obtained  easements  to 
prevent  the  draining  or  filling  in  of  wetland  and  that  included 
the  surrounding  private  land  in  the  tract;  and  private  land 
that  included  private  land  and  other  state  and  federal  land 
not  owned  by  the  service.  Townships  were  assigned  to  3 
strata  by  the  following  rules: 

(1)  A  low  landownership  stratum  contained  15.5  km  or 
less  of  easements  and  0  km^  of  land  owned  by  the 
service. 

(2)  A  moderate  landownership  stratum  contained  greater 
than  15.5  km^  of  easements  and  less  than  2.6  km^  of 
land  owned  by  the  service. 

(3)  A  high  landownership  stratum  contained  2.6  km^  or 
more  of  land  owned  by  the  service. 

We  randomly  drew  10.4-km^  plots  from  each  land- 
ownership  stratum  to  obtain  a  sample  with  higher  sampling 
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Fig.  1.  Distribution  of  10.4-kni^  sample 
plots  ■  used  to  evaluate  duck 
(Anatinae)  habitat  and  population 
size  during  1987-90  in  the  prairie 
pothole  region  of  the  United  States. 
Waterfowl  management  districts 
used  for  analysis  of  duck  popula¬ 
tions  are  shaded.  Letter  codes  are: 
AR,  Arrowwood;  AU,  Audubon; 
CL  Crosby-Lostwood;  DL,  Devils 
Lake;  JC,  J.  Clark  Salyer;  KU, 
Kulm;  LL,  Long  Lake;  ML,  Medi¬ 
cine  Lake;  TE,  Tewaukon;  WB, 
Waubay. 


fraction  in  areas  with  high  service  landownership  or  ease¬ 
ment  because  these  areas  were  most  desirable  for  simula¬ 
tions  of  management.  At  that  time,  recent  color-infrared 
aerial  photographs  of  only  422  of  the  selected  500  plots  were 
available.  The  photographs  were  taken  in  May  of  a  wet  year. 
The  final  sample  had  sampling  fractions  of  0.0045  in  the 
low,  0.01 00  in  the  moderate,  and  0.0824  in  the  high  service- 
landownership  strata.  Although  considered  sufficiently  rep¬ 
resentative  for  model  simulations,  the  sample  was  no  longer 
strictly  random  because  of  the  plots  without  photographs 
(Cowardin  et  al.  1988). 

We  used  the  existing  422  sample  plots  as  the  basis  for  the 
remote-sensing-based  system  because  of  the  large  prior 
expenditure  for  mapping  and  digitizing  data  from  those 
plots.  This  decision  created  three  problems.  First,  the  re¬ 
mote-sensing-based  system  required  estimates  by  wetland 
management  districts  that  were  not  considered  in  the  original 
sample  selection.  Second,  the  method  of  stratification  re¬ 
sulted  in  some  plots  without  service  landownership  in  the 
high  or  moderate  landownership  strata  because  townships 


were  assigned  to  strata  and  the  rules  did  not  apply  to  the 
plots.  Third,  when  the  wetland-management-district 
boundaries  were  placed  over  the  existing  sample,  some 
waterfowl  management  districts  contained  few  sample  plots 
and  some  strata  inside  waterfowl  management  districts  were 
not  represented. 

We  restratified  the  sample  to  overcome  the  first  two 
problems  by  estimating  landownership  in  each  10.4-km“ 
cell  of  the  sample  universe  and  by  assigning  each  plot  to  a 
wetland  management  district.  This  procedure  allowed  us  to 
calculate  weights  for  each  stratum  and  thus  obtain  unbiased 
estimates  of  each  parameter  in  each  waterfowl  management 
district.  The  restratification  was  accomplished  by  mapping 
and  digitizing  all  landownership  classes  in  the  entire  sample 
universe  on  1 :250,000  USGS  maps  and  by  then  overlaying 
a  grid  of  all  10.4-km^  plots.  Where  wetland-management- 
district  boundaries  lay  along  rivers,  the  areas  were  divided 
into  irregularly  shaped  plots  of  approximately  10.4  km  . 
This  grid  was  digitized,  and  the  digital  data  were  merged 
with  the  landownership  map  by  The  Map  Overlay  and 


Table  1.  Numbers  of  10.4-km^  plots  (n)  for  evaluating  duck  (Anatinae)  habitat  and  estimating  numbers  of  ducks  and 
number  of  plots  covered  by  videography  (NV)^  during  1987-90  in  the  prairie  pothole  region  of  the  United  States. 


State 

1987 

1988 

1989 

1990 

N 

NV 

N 

NV 

N 

NV 

N 

NV 

Minnesota 

95 

87 

98 

79 

128 

118 

128 

128 

Montana 

14 

14 

14 

13 

14 

14 

14 

14 

North  Dakota 

203 

202 

226 

219 

226 

220 

226 

223 

South  Dakota 

23 

23 

23 

22 

75 

74 

75 

75 

All  states 

335 

326 

361 

333 

443 

426 

443 

440 

^Weather  conditions  and  time  sometimes  prevented  aquiring  video  data  in  the  required  time  interval. 
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Fig.  2.  Regression  of  number  of  blue¬ 
winged  teal  {Anas  discors),  gad- 
wall  (Anas  strepera)^  mallard 
(Anas  platyrhynchos),  northern 
shoveler  (Anas  clypeata),  and 
northern  pintail  (Anas  acuta)  pairs 
on  pond  size  during  1987-90  in 
the  prairie  pothole  region  of  the 
United  States. 


Statistical  System  (Pywell  and  Niedzweadk  1 980)  to  pro¬ 
duce  a  file  with  the  landownership  of  each  1 0.4-km^  cell  in 
the  universe.  Based  on  these  data,  all  plots  in  the  universe 
were  assigned  to  strata  by  the  following  rules: 

1.  A  refuge  stratum  included  any  plot  that  contained  any 
national  wildlife  refuge  land  regardless  of  other  land- 
ownership  in  the  plot. 

2.  A  waterfowl-production-area  stratum  included  any  plot 

not  classified  as  a  refuge  that  included  64.8  ha  or  more 
of  waterfowl  production  area. 

3.  An  easement  stratum  included  any  plot  not  classified  as 

refuge  or  waterfowl  production  area  and  64.8  ha  or 
more  of  service  easement. 

4.  A  private  stratum  included  any  plot  not  in  the  previous 
three  strata. 

In  addition,  each  plot  in  the  universe  was  assigned  to  the 
appropriate  wetland  management  district. 

To  overcome  the  third  problem,  we  added  additional 
plots  inside  waterfowl  management  districts  that  had  insuf¬ 
ficient  sample  size.  We  required  at  least  two  sample  plots  in 
each  landownership  in  each  wetland  management  district. 
The  sample  contained  few  plots  from  the  refuge  stratum. 
Refuges  are  often  larger  than  10.4-km^  plots,  and  the  plot 
size  is  not  well  suited  to  obtaining  a  representative  sample. 
Therefore,  for  this  report,  we  collapsed  refuge  (stratum  1) 
and  waterfowl  production  area  (stratum  2)  into  a  single 
stratum  called  service. 

Sample  Wetland  Basins 

Approximately  200  wetland  basins  were  selected  from 
all  plots  in  each  wetland  management  district  as  sites  for 


conducting  breeding-pair  counts.  Sample  size  was  deter¬ 
mined  according  to  the  amount  of  time  available  for  con¬ 
ducting  pair  counts  rather  than  by  statistical  estimation  of 
sample  size  required  for  a  given  degree  of  precision.  The 
purpose  of  the  pair  counts  was  to  adjust  baseline  regressions 
for  annual  and  geographic  variation  in  pair  density.  There¬ 
fore,  we  used  an  optimum  allocation  for  a  stratified  random 
sample  and  treated  the  wetland-basin  classes  as  strata  to 
obtain  a  sample  throughout  the  range  of  wetland-basin  sizes 
and  to  avoid  oversampling  of  small  basins  that  are  often  dry 
and,  therefore,  provide  no  information  about  duck  density. 
Although  the  technique  is  intended  for  minimizing  the  vari¬ 
ance  of  a  stratified  mean,  given  the  individual  strata  vari¬ 
ances  (see  Neyman  allocation  in  Cochran  1977),  it  also  had 
the  desired  effect  of  reducing  the  sample  of  temporary 
wetlands  basins,  which  had  a  smaller  variance  than  the  more 
permanent  wetland-basin  classes.  Strata  variances  were  es¬ 
timated  from  a  regression  model  by  obtaining  estimates  of 
the  number  of  mallard  pairs  in  each  wetland  basin  in  each 
plot.  The  area  of  each  wetland  basin  was  obtained  from 
special  maps  prepared  by  the  National  Wetland  Inventory. 
All  wetland  basins  were  assumed  to  contain  ponds.  Vari¬ 
ances  of  the  number  of  mallard  pairs  among  wetland  basins 
within  each  wetland-basin  class  in  each  wetland  manage¬ 
ment  district  were  then  calculated. 

Habitat  Data 

Classifications 

We  required  a  classification  of  wetland  basins  to  estimate 
breeding  duck  populations  and  a  classification  of  upland 
and  wetland  nesting  habitat  to  estimate  duck  production. 
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Wetland  was  treated  different  from  the  other  habitat  classes. 
The  National  Welland  Inventory  mapped  wetland  and  up¬ 
land  on  the  plots  as  a  special  project.  Welland  was  mapped 
according  to  the  classification  and  definitions  in  Cowardin 
et  al.  (1979)  and  with  some  exceptions  according  to  the 
current  mapping  conventions  by  the  National  Wetland  In¬ 
ventory.  The  exceptions  were  that  no  codes  for  unknown 
water  regime  or  mixed  classes  were  allowed.  Wetland  map¬ 
ping  of  the  plots,  except  the  addition  of  a  unique  number  for 
each  polygon  in  each  basin,  was  essentially  identical  to  the 
National  Wetland  Inventory  operational  inventory  (Wilen 
1990).  Cowardin  (1982)  illustrated  the  difference  between 
classifications  of  wetland  area  and  of  wetland  basins.  Avail¬ 
able  data  for  constructing  pair-wetland  regressions  were 
from  various  basin  classifications.  Therefore,  we  had  to 
translate  the  cover  classes  mapped  by  the  National  Wetland 
Inventory  into  basin  classes.  The  technique  was  a  simplifi¬ 
cation  of  the  rules  by  Stewart  and  Kantrud  ( 1 97 1 )  for  form¬ 
ing  pond  classes  from  wetland  zones.  Their  pond  classes 
(equivalent  to  our  basin  classes;  Table  2)  were  derived  from 
the  water  regime  of  the  zone  with  most  permanence  and  with 
an  aerial  cover  of  5%  or  more.  Our  algorithm  first  summed 
the  area  of  all  wetland  polygons  in  a  basin  by  a  unique 
identifier  coded  at  the  time  of  digitization.  It  then  searched 
for  the  polygon  with  the  most  permanent  water  regime 
regardless  of  polygon  size.  If  two  or  more  polygons  had  the 
same  water  regime,  the  algorithm  selected  the  largest.  That 
polygon  became  the  basin-naming  polygon  and  was  used  to 
determine  the  class  of  the  basin  (Table  2).  If  the  basin 
contained  only  one  wetland  polygon,  that  polygon  became 
the  basin-naming  polygon  by  default. 

Upland  habitats  were  those  described  by  Cowardin  et  al. 
(1988):  grassland,  hayland,  planted  cover,  cropland,  wood¬ 
land,  scrubland,  other,  right-of-way,  and  barren  land.  We 
modified  some  of  the  mapped  habitats  to  obtain  the  nesting- 
habitat  classes  based  on  landownership  and  to  reflect  major 
habitat  changes  since  mapping.  Grassland  was  subdivided 
into  grassland  and  grassland-wildlife.  Grassland  was  typi- 


Table  2.  Wetland-basin  class  definitions  based  on  the 
water  regime  of  the  basin-naming  polygon.^ 


Basin  class 

Basin-Naming  Polygon 

Temporary  basin 

Water  regime  temporary  (a)*^  or 
saturated  (b) 

Seasonal  basin 

Water  regime  seasonal  (c) 

Semipermanent  basin 

Water  regime  semipermanent  (0 

Lake 

System  Lacustrine  (L)  or  water 
regime  intermittently  exposed  (h) 
or  permanent  (g) 

^The  basin-naming  polygon  is  the  polygon  with  the  most  permanent 
water  regime  in  a  wetland  basin. 

^Letters  in  parentheses  refer  to  symbols  on  National  Wetland  Inventory 
maps. 


cal  pasture.  Grassland-wildlife  was  nonuse  grassland  that  is 
typical  in  waterfowl  production  areas  and  on  refuges.  Be¬ 
cause  these  types  are  difficult  to  separate  on  aerial  photo¬ 
graphs,  we  defined  any  grassland  in  waterfowl  production 
areas  or  refuges  as  grassland-wildlife.  The  Conservation 
Reserve  Program  of  the  Food  Security  Act  of  1985  (United 
States  Congress  1985)  went  into  effect  during  our  develop¬ 
ment  of  the  remote-sensing-based  system.  Locations  of 
Conservation  Reserve  Program  tracts  on  10.4-km^  plots 
were  furnished  by  wetland-managcment-district  managers 
who  obtained  the  locations  from  the  county  offices  of  the 
Agricultural  Stabilization  and  Conservation  Service  in  each 
wetland  management  district.  We  used  the  class  Consen^a- 
tion  Reserve  Program  cover  for  land  in  these  tracts.  Prior  to 
analysis,  we  modified  the  corresponding  habitat  files  for 
wetland  management  districts  where  we  estimated  duck 
production.  The  amounts  of  planted  cover  in  waterfowl 
production  areas  changed  from  the  time  photographs  for  the 
National  Wetland  Inventory  were  taken  and  the  time  the 
plots  were  mapped.  We  obtained  maps  of  current  planted 
cover  tracts  in  each  wetland  management  district  from  the 
managers  in  1987  and  modified  our  data  files  to  reflect 
current  conditions.  We  used  these  data  for  all  years  irrespec¬ 
tive  of  possible  minor  changes  during  1987-90. 

Habitat  Maps 

Habitat  on  each  plot  was  interpreted,  mapped,  and  digit¬ 
ized  by  the  National  Wetland  Inventory.  The  mapping  was 
a  special  project  conducted  for  us  prior  to  production  of 
standard  wetland  maps  by  the  National  Wetland  Inventory. 
Data  were  from  high-altitude  ( 1 :63,360),  color-infrared  pho¬ 
tographs  taken  during  the  late  1970’s  and  early  1980’s.  All 
features  on  the  plots  were  delineated  with  a  5-aught  pen  on 
acetate  overlays.  Wherever  possible,  areas  were  shown  as 
closed  polygons,  but  some  features  such  as  roads,  trails,  and 
rock  piles  had  to  be  shown  as  lines  or  points  because  of  the 
small  scale.  At  the  time  of  digitization,  a  unique  basin 
number  was  added  as  an  attribute  to  all  polygons  in  a  single 
wetland  basin.  Polygons,  linear  features,  and  points  were 
transferred  to  1 :24,000  USGS  quadrangle  maps  by  a  Bausch 
and  Lomb  Zoom  Transfer  Scope. 

The  resulting  maps  were  then  digitized  on  a  digitizing 
tablet  and  converted  to  Map  Overlay  and  Statistical  System 
files  (Pywell  and  Niedzweadk  1980).  A  second  set  of 
1 :24,000  plot  maps  showing  landownership  boundaries  was 
prepared  from  data  on  file  at  realty  offices  of  the  U.S.  Fish 
and  Wildlife  Service.  These  maps  were  also  digitized  into 
Map  Overlay  and  Statistical  System  files.  The  two  Map 
Overlay  and  Statistical  System  files  were  overlaid  to  pro¬ 
duce  a  file  with  the  landownership  attributes  of  all  polygons. 
Text  files  of  each  polygon  with  single  records  were  produced 
from  these  files. 
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Because  the  remote-sensing-based  system  required  that 
all  habitats  have  some  area,  line  and  point  features  were 
buffered  by  multiplying  length  by  average  feature  width 
(linear  features)  or  by  nt^  (point  features).  The  following 
average  dimensions  were  determined  from  aerial  photo¬ 
graphs  and  were  used  for  buffering  linear  and  point  data 
on  maps:  8.2  m  wide  for  shelterbelts;  14.6  m  diameter  for 
rockpiles  and  brush  or  grass  areas,  too  small  to  enclose 
with  a  polygon;  1 4.6  m  wide  for  linear  wetland  basins;  and 
15.3  m  diameter  for  point  wetland  basins.  Linear  road 
features  were  buffered  for  the  width  of  the  road  surface, 
which  we  equated  to  being  barren,  and  for  the  distance 
from  the  road  surface  to  the  border  of  the  right-of-way. 
Distances  from  the  center  line  of  the  road  were  3.1  m  to 
the  edge  of  the  road  surface  and  10.1  m  to  the  far  edge  of 
the  right-of-way  on  gravel  roads  and  respectively  3.8  m 
and  19.1  m  on  hard-surface  roads  and  6.1  m  and  19.8  m 
on  railroads.  The  average  width  of  fence  rows  and  vege¬ 
tated  strips  between  fields  was  3.1  m. 

Areas  for  buffered  linear  features  and  points  were  added 
to  the  text  files  derived  from  polygon  data.  This  inflated  the 
total  plot  area.  All  polygon  areas  were  then  rescaled  to  the 
true  plot  area  by  calculating  a  correction  factor  (true  plot 
area  -5-  inflated  plot  area)  and  by  multiplying  the  area  of  each 
polygon  by  the  correction  factor. 

Aerial  Videography 

The  maps  and  databases  contained  data  on  all  wetland 
basins  (wet  or  dry)  when  the  photographs  were  taken. 
We  assumed  that  the  maps  had  no  errors  from  omissions 
or  commissions.  The  remote-sensing-based  system  also 
required  that  we  know  the  numbers  and  sizes  of  all  ponds 
(wetland  basins  with  water)  each  spring.  We  selected 
aerial  video  taken  during  flights  in  early  May  of  each 
year  as  the  technique  to  obtain  these  data.  Video,  al¬ 
though  lacking  the  good  resolution  of  photographs,  has 
advantages  over  photographs  (Sidle  and  Ziewitz  1990). 
Video  is  less  expensive  than  aerial  photography.  Be¬ 
cause  a  monitor  is  in  the  aircraft,  the  user  can  guide  the 
pilot  to  the  target  area  and  knows  whether  the  target  area 
was  recorded.  The  data  are  ready  for  immediate  use  at 
the  completion  of  the  flight.  The  Map  and  Image  Proc¬ 
essing  System  software  (Miller  et  al.  1 990)  allows  instant 
capture  of  the  images  in  digital  form  directly  from  the 
video  signal  (unlike  photographs,  which  must  be 
scanned  to  produce  digital  data). 

We  used  a  Panasonic  D  5000  video  camera  with  a 
5.9-mm  Angenieux  lens,  a  Panasonic  AG  2400  video  re¬ 
corder,  and  a  Panasonic  CT  500V  monitor.  The  camera  was 
controlled  by  a  Panasonic  WVCR12  controller  board.  The 
camera  was  mounted  in  a  custom-designed  aluminum 
mount  that  allowed  the  leveling  and  rotation  of  the  camera 
to  correct  for  crabbing  of  the  aircraft.  The  use  of  the  short 


lens  did  not  allow  use  of  the  automatic  white-balance  and 
iris.  The  aperture  had  to  be  set  before  the  flight,  and  the 
camera  was  white-balanced  just  before  or  during  the  flight. 

We  used  several  aircraft  including  Cessnas  172, 
172RG,  and  182  and  a  Partenavia  Surveyor.  A  variety  of 
camera  port  sizes  and  locations  was  used.  Camera  ports 
from  12.7  to  30.5  cm  in  diameter  proved  adequate,  but 
larger  ports  were  easier  to  use  because  the  chance  of 
including  part  of  the  aircraft  skin  in  the  video  frame  is 
less.  Video  data  were  obtained  from  an  altitude  of 
3,812.5  m  above  ground  level.  This  elevation  combined 
with  our  type  of  lens  permitted  a  5.2-km  swath  width  and 
adequate  room  for  navigational  error. 

Analysis  of  Image  Data 

Video  data  were  captured  in  digital  form  by  a  microcom¬ 
puter  with  aTarga  1 6  image-capture  board  and  The  Map  and 
Image  Processing  System  (MIPS)  software.  This  procedure 
produced  a  16-bit  composite  raster.  After  capture,  we  clas¬ 
sified  ponds  by  the  Feature  Mapping  procedure  in  MIPS. 
The  objective  was  to  classify  all  areas  that  are  covered  by 
water,  including  vegetation  growing  in  water.  The  proce¬ 
dure  required  skill  in  interpretation  and  knowledge  of  local 
wetland  conditions.  The  percent  of  the  wetland  basin  cov¬ 
ered  by  water  was  recorded  during  the  counting  of  the 
breeding  pairs  and  furnished  ground  truth.  Feature  Mapping 
in  MIPS  can  be  used  in  either  an  automated  mode  or  by 
drawing  boundaries  of  a  pond  on  the  screen  with  a  mouse. 
Video  data  seldom  furnish  sufficient  spectoral  separation  for 
completely  automated  classification  of  a  scene.  We  picked 
and  classified  training  pixels,  known  to  contain  water  from 
ground  observation,  until  errors  of  commission  began  to 
appear.  It  was  then  necessary  to  begin  on-screen  interpreta¬ 
tion  by  drawing  boundaries  around  areas  interpreted  as  wet. 
Where  the  basin  contained  emergent  vegetation,  we  looked 
for  water  along  the  shore  or  in  openings  in  the  vegetation. 
This  was  a  good  indication  of  vegetation  underlain  by  water. 
Sun-glint  problems  were  resolved  by  referring  to  the  original 
video  tape  and  by  observing  sun-glint  move  across  the  scene 
as  the  aircraft  moved  over  the  wetland  basin.  Cloud-shadow 
problems  were  overcome  by  comparing  the  relation  between 
the  darker  shading  of  water  and  the  lighter  shading  of  upland 
in  shadows  and  clear  areas. 

Interpretation  is  subject  to  errors,  and  consistency  among 
interpreters  is  important.  Two  people  interpreted  the  same 
scenes  of  most  video  data  to  identify  errors  in  the  classifica¬ 
tion  of  amount  of  water.  When  inconsistencies  occurred,  the 
area  was  reclassified. 

Duck  Populations 

Study  Area 

We  restricted  our  analysis  of  annual  change  in  ponds, 
duck  population  sizes,  and  duck  production  to  10  waterfowl 
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management  districts  (Fig.  1 )  for  which  we  had  aerial  video 
data  during  1 987-90.  Although  video  was  obtained  of  other 
waterfowl  management  districts  as  they  were  added  to  the 
remote-sensing-based  system,  the  lack  of  data  in  some  years 
would  have  confounded  comparisons  among  years. 

Regression  Estimates 

Estimates  of  breeding-population  sizes  were  derived 
from  a  regression  model.  Data  for  constructing  baseline 
regressions  for  all  species  except  wood  ducks  (Aix 
sponsa)  were  gathered  at  the  Arrowwood  Waterfowl 
Management  District  (1982-84).  For  wood  ducks,  we 
used  data  gathered  at  the  Fergus  Falls  (1986-87)  and 
Detroit  Lakes  (1987)  waterfowl  management  districts. 
Baseline  regressions  were  constructed  by  regressing  the 
number  of  pairs  observed  on  each  pond  on  the  size  and 
the  square  root  of  the  size  of  the  pond.  The  general  form 
of  the  equations  was: 

Pairs=Axwet  area  +  Bx-^wet  area  (1) 

Wet  area  in  this  equation  is  the  wet  area  in  each  pond. 
(Estimates  of  the  regression  coefficients  A  and  B  are 
presented  in  Table  3,  and  the  regression  curves  are  pre¬ 
sented  on  Fig.  2). 

Ground  Counts 

The  boundary  of  each  wetland  basin  where  breeding-pair 
counts  were  conducted  was  delineated  on  aerial  photo¬ 
graphs,  and  a  unique  identification  number  for  the  naming 
of  the  polygon  of  that  basin  was  placed  on  the  photograph. 

Table  3.  Baseline  regression  coefficients  used  for 
estimating  breeding  duck  (Anatinae)  populations  on 
10.4-km^  plots  from  the  area  (ha)  of  individual  ponds 
determined  from  airborne  video  in  the  prairie  pothole 
region  of  the  United  States,  1987-90. 

Regression  coefficients^ 

Species  A  B 


Mallard 

0.0106 

0.2899 

Gadwall 

0.0341 

0.2848 

American  wigeon 

0.0000 

0.0193 

Green-winged  teal 

0.0000 

0.0431 

Bliie-winged  teal 

0.0000 

0.7376 

Northern  shoveler 

0.0136 

0.1870 

Northern  pintail 

0.0000 

0.1866 

Redhead 

0.0410 

0.2233 

Canvasback 

0.0064 

0.0453 

Lesser  scaup 

0.0173 

0.0528 

Ring-necked  duck 

0.0022 

0.0110 

Ruddy  duck 

0.0119 

0.1226 

Wood  duck 

0.0064 

0.5612 

^  Pairs  =  A  X  wet  area  +  Bx  wet  area. 


The  pair  count  was  conducted  by  a  person  who  walked 
around  each  pond.  Pairs  on  lakes  and  rivers  were  sometimes 
surveyed  from  a  boat  or  from  a  canoe.  When  dense  emergent 
vegetation  was  present,  the  observer  moved  through  the 
vegetation  and  flushed  the  birds.  Sample  ponds  were  more 
widely  dispersed  in  the  remote-sensing-based  system  than 
in  surveys  that  are  often  conducted  on  plots  or  transects.  We 
assumed  that  sample  ponds  were  usually  far  enough  apart  to 
prevent  the  flushing  of  birds  from  one  sample  pond  to 
another.  Social  groups  were  recorded  on  field  forms  de¬ 
signed  for  data  entry  into  a  microcomputer.  Data  were 
entered  soon  after  completion  of  the  counts.  Error  checking 
was  done  by  a  series  of  custom-designed  programs  that  were 
executed  at  the  time  of  data  entry  and  immediately  after  data 
entry.  The  size  of  a  breeding  population  was  expressed  as 
the  estimated  number  of  breeding  pairs,  which  included 
observed  breeding  pairs  and  pairs  that  were  inferred  from 
other  social  groups.  We  followed  approximately  the  same 
techniques  described  by  Hammond  (1969)  and  Dzubin 
(1969).  The  estimated  number  of  breeding  pairs  was  esti¬ 
mated  by  the  computer  at  the  time  of  data  entry  as  follows: 
1  observed  pair  =  1  estimated  breeding  pair;  1  lone  male  = 
1  estimated  breeding  pair;  1  lone  female  diving  duck  =  1 
estimated  breeding  pair;  each  flocked,  male  dabbling  duck 
in  flocks  of  five  or  fewer  individuals  =  1  estimated  breeding 
pair  except  in  wigeons  and  in  northern  shovelers  for  which 
the  number  of  flocked  males  was  not  converted  to  estimated 
number  of  breeding  pairs;  each  flocked  female  diving  duck 
in  flocks  of  five  or  fewer  individuals  =  1  estimated  breeding 
pair.  Flocked  birds  in  flocks  with  more  than  five  individuals 
were  not  used  to  estimate  the  number  of  breeding  pairs. 

Data  Analyses 

Estimation  Procedures 

We  required  two  types  of  estimates:  totals  (e.g.,  num¬ 
ber  of  ponds)  and  ratios  of  totals  (e.g.,  average  pond  size 
=  area  of  ponds  per  number  of  ponds).  Where  possible, 
we  also  wanted  variances  of  these  estimates. 

2 

To  estimate  totals,  we  treated  our  sample  of  10.4-km 
plots  as  a  random  sample  and  stratified  by  waterfowl  man¬ 
agement  district  and  U.S.  Fish  and  Wildlife  Service  land- 
ownership.  The  amount  of  area  in  a  given  landownership 
varied  from  plot  to  plot.  We  employed  a  combined  ratio 
estimator  (Cochran  1977)  that  made  use  of  the  fact  that  the 
totals  that  we  wanted  to  estimate  positively  correlated  with 
the  amount  of  area  in  that  landownership.  Estimation  was 
the  computation  of  an  average  density  by  landownership 
from  the  sample  plots  in  each  waterfowl  management  dis¬ 
trict  and  the  multiplication  of  that  value  by  the  total  amount 
of  area  of  that  landownership  in  each  waterfowl  manage¬ 
ment  district.  The  overlay  of  plots  on  1:250,000  land- 
ownership  maps,  described  under  stratification,  provided 
estimates  of  the  amount  of  land  by  landownership  in  each 
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stratum  and  in  each  waterfowl  management  district.  The 
combined  totals  of  all  landownerships  were  estimated  by 
summing  the  landownership  totals.  Totals  by  state  and  their 
estimated  variances  were  calculated  by  summing  estimates 
from  constituent  waterfowl  management  districts. 

Variances  of  the  landownership  totals  were  calculated 
with  the  Formula  6.51  from  Cochran  (1977).  An  upper 
bound  to  the  variance  of  the  combined  total  landownership 
of  all  ownership  was  estimated  from  Formula  5.10  (Cochran 
1977:93). 

We  divided  appropriate  estimates  of  population  totals  of 
one  parameter  by  totals  of  another  parameter  to  arrive  at 
estimates  of  population  ratios  (e.g.,  pairs  per  pond).  We  do 
not  report  estimated  variances  of  population  ratios  because 
exact  formulas  do  not  exist  and  formulas  for  approximations 
are  complex  (Mood  et  al.  1974). 

Estimates  of  Habitat  Parameters 

Number  of  wetland  basins,  area  of  wetland  habitat, 
and  amount  of  upland  nesting  cover  were  obtained  from 
the  maps  of  10.4-km^  plots  prepared  by  the  National 
Wetland  Inventory.  The  number  of  ponds  and  the  area  of 
wetland  covered  by  water  varied  from  year  to  year  and 
were  obtained  from  the  classified  video  scenes. 

Number  of  Breeding  Pairs 

To  estimate  the  size  of  the  breeding  duck  population,  the 
area  of  each  pond  (from  the  video  scenes)  was  entered  into 
the  baseline  regression  of  each  species  to  estimate  the  num¬ 
ber  of  the  breeding  pairs.  These  by-pond  estimates  were 
summed  by  plot  and  species  and  then  expanded  to  waterfowl 
management  districts  as  explained  undox  Estimation  Proce¬ 
dures.  The  resulting  estimates  were  based  on  the  assumption 
that  the  densities  of  ducks  on  ponds  of  the  same  size  remains 
constant  from  year  to  year  and  from  area  to  area.  An  adjust¬ 
ment  of  the  differences  among  areas  and  years  in  each 
waterfowl  management  district  was  based  on  counts  of  the 
sample  wetland  basins.  A  correction  (y)  was  defined  as: 

_ _ Total  number  of  counted  pairs _ 

^  Number  of  pairs  predicted  by  regression  ‘ 

(2) 

To  determine  the  number  of  predicted  pairs  in  this 
equation,  we  required  the  area  of  each  counted  pond.  The 
video  data  did  not  identify  each  pond  by  number. 
Therefore,  we  used  the  product  of  the  estimated 
proportion  of  the  basin  that  contained  water  and  the 
mapped  area  of  the  basin  to  estimate  the  area  of  each 
counted  pond.  Finally,  the  estimated  number  of  predicted 
pairs  in  the  waterfowl  management  district  was 
multiplied  by  yto  correct  for  differences  in  pair  densities 
among  years  and  among  waterfowl  management 
districts. 


Estimation  of  Production 

The  number  of  young  ducks  recruited  to  the  population 
in  fall  cannot  be  counted  and,  therefore,  we  used  a  model  to 
calculate  estimates.  The  model  was  deterministic  and  pro¬ 
duced  a  result  similar  to  that  obtained  from  a  stochastic 
mallard  productivity  model  (Johnson  et  al.  1987).  For  our 
purposes,  the  deterministic  technique  was  simpler  and  per¬ 
mitted  more  rapid  calculation  of  estimates  than  the  stochas¬ 
tic  model.  An  example  illustrates  the  method.  To  estimate 
the  number  of  recruits  in  each  landownership  class,  we 
essentially  worked  backwards  from  the  general  relation: 

Recruits  =  2Rn,  (3) 

where  2  is  a  constant  based  on  the  assumption  of  equal  sex 
ratio  at  hatch,  R  is  the  recruitment  rate  as  defined  by 
Cowardin  and  Johnson  (1979),  and  n  is  the  number  of 
breeding  pairs. 

The  recruitment  rate  was  broken  down  into  the  compo¬ 
nents  described  by  Cowardin  and  Johnson  (1979)  in  the 
equation: 


where  H  is  hen  success,  2  is  a  constant  used  because  recruit¬ 
ment  is  defined  in  terms  of  only  females,  Z  is  the  proportion 
of  entire  broods  that  survived  to  time  of  census,  and  B  is  the 
average  brood  size  at  fledging.  For  mallards,  we  used  Z  = 
0.74  and  B  =4.9  from  Cowardin  and  Johnson  (1979).  Brood 
survival  (Z)  of  0.74  was  also  used  for  gadwalls,  blue- winged 
teals,  northern  shovelers,  and  pintails.  The  mallard  brood 
size  of  4.9,  adjusted  for  broods  of  size  zero  (Cowardin  and 
Johnson  1979),  was  used  as  baseline  to  estimate  Class-m 
brood  (Gallop  and  Marshal  1954)  size  of  other  species 
(Table  4). 

Hen  success  was  determined  from  the  relation: 

=  (5) 

where  P  is  clutch  success,  a  is  an  index  to  nesting  effort 
(Cowardin  and  Johnson  1979),  and  e  is  approximately 
2.718.  A  relation  between  a  and  the  percentage  of  basins 
containing  water  was  derived  from  unpublished  data  gath¬ 
ered  during  a  study  of  mallards  in  central  North  Dakota 
(Table  5).  The  technique  was  the  same  as  that  used  by 
Cowardin  et  al.  (1985)  for  relating  a  to  the  percentage  of 
semipermanent  ponds  containing  water.  The  resulting  equa¬ 
tion  was: 

a=  0.616  +  0.00603  (W),  (6) 

where  W  is  the  percentage  of  basins  that  contained  water 
in  a  waterfowl  management  district,  estimated  from  the 
video  data  and  divided  by  total  number  of  basins  as 
mapped  by  the  National  Wetland  Inventory. 
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Table  4.  Average  brood  sizes  of  surface-feeding  ducks  (Anatinae),  derived  by  scaling  mallard  {Anas  piatyrhynchos) 
brood  size  by  relative  clutch  sizes  of  mallards  and  other  species  during  a  1 987-90  study  in  the  prairie  pothole  region 
of  the  United  States. 


Species 

Clutch  size^ 

Average  brood  size^ 

Mallard 

8.4 

4.90 

Gadwall 

9.9 

5.78 

Blue-winged  teal 

10.2 

5.93 

Northern  shoveler 

9.9 

5.78 

Northern  pintail 

7.1 

4.12 

^ After  Duebbert  and  Frank  (1 984: Table  3) 

Average  brood  size 

=  clutch  size  ^  mallard  clutch  size  x  4.90 

Table  5.  Percent  of  total  wetland  basins  that  contained  ponds  in  a  study  area  in  North  Dakota  (Cowardin  et  al.  1985) 

and  estimates  of  during  1977- 

-80  in  the  prairie  pothole  region  of  the  United  States. 

Year 

Total  basins 

Ponds 

% 

a 

1977 

1984 

224 

11.3 

0.5918 

1978 

1984 

880 

44.4 

0.9108 

1979 

1984 

1527 

77.0 

1,0444 

1980 

1984 

507 

25.6 

0.8700 

®  An  index  to  nesting  intensity  after  Cowardin  and  Johnson  (1979). 


The  number  of  produced  recruits  on  a  plot  can  be  deter¬ 
mined  from  the  previous  two  equations  if  clutch  success  in 
each  plot  is  available.  Clutch  success  and  the  allocation  of 
nests  to  the  landownership  classes  were  computed  as  in  the 
following  simplified  example. 

We  assumed  that  the  distribution  of  habitat  and  land- 
ownership  (Fig.  3)  was  two  habitats  A  and  B  and  two 
landownerships,  259.2  ha  of  service  landownership  and 


777.6  ha  of  private  landownership.  We  assumed  that  the 
ducks’  preference  for  habitat  A  was  twice  that  for  habitat  B. 
If  n  nests  were  on  64.8  ha  of  habitat  B,  the  number  of  nests 
on  the  other  tracts  of  land  could  be  calculated  as  area  of  the 
tract  64.8  X  preference  value  (Fig.  3).  If  clutch  success 
was  assumed  to  be  Pa  =  0.50  and  Pb  =  0.10,  we  could 
calculate  the  number  of  successful  nesting  attempts  in  habi¬ 
tats  A  and  B  (Fig.  4). 


259.2  ha 

H 

Habitat  =  A 

8n  nests 

518.4  ha 

Habitat  =  B 
8n  Nests 


Fig.  3.  Hypothetical  distribution  of  duck  (Anat- 
inae)  habitats  and  nests  on  a  10.4-km^  plot 
with  two  landownerships  and  two  habitats. 
Landownership  classes  include  land  owned 
by  the  U.S.  Fish  and  Wildlife  Service  (serv¬ 
ice)  and  privately  owned  land  (private)  in  the 
prairie  pothole  region  of  the  United  States, 
1987-90.  A  and  B  denote  arbitrary  land 
cover  types. 


/ 


Fws  m 


Private  O 
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Fig.  4.  Numbers  of  successful  nesting  attempts 
by  ducks  (Anatinae)  on  a  hypothetical  10.4- 
km^  plot  with  two  landownerships  and  two 
habitats.  Landownership  classes  include  land 
owned  by  the  U.S.  Fish  and  Wildlife  Service 
(service)  and  privately  owned  land  (private) 
in  the  prairie  pothole  region  of  the  United 
States,  1987-90. 


With  these  assumptions,  the  total  number  of  nests  (Fig.  3) 
is  23/7  and  the  number  of  successful  nesting  attempts  (Fig.  4) 
is  7.9/7.  Therefore,  the  weighted  clutch  success  in  the  entire 
plot  (P)  is: 

^  N umber  of  successful  nesting  attempts 

~  Total  nesting  attempts 

=  2^=0. 34. 

23n  (7) 

To  allocate  recruits  to  the  landownership  classes,  we 
calculated  the  proportion  of  clutches  hatched  in  each  land- 
ownership  class. 

Results  and  Discussion 

Habitat  Base 

The  amount  of  land  in  the  sample  universe  was 
3 15,980  km^.  Most  land  was  in  North  Dakota  (42%),  and 
more  land  was  in  South  Dakota  (29%)  than  in  Minnesota 
(25%)  and  in  Montana  (5%).  These  values  reflected  the 


proportions  of  the  prairie  pothole  region  in  each  state 
except  in  Montana  where  our  data  were  taken  from  only 
one  waterfowl  management  district  in  the  northeastern 
part  of  the  state  and  in  Iowa  a  state  that  was  not  included 
in  the  survey.  Most  of  the  land  in  each  state  was  in  private 
landownership  (Table  6).  The  difference  among  states 
reflected  differences  in  the  acquisition  of  easements 
rather  than  variation  in  the  amount  of  wetland.  Only 
about  1%  of  the  land  in  the  prairie  pothole  region  is  in 
service  landownership.  This  land,  including  wetland  and 
upland,  has  the  greatest  potential  for  intentional  modifi¬ 
cation  of  habitat  to  benefit  ducks  because  the  service 
controls  management. 

Wetland  Habitat 

We  estimated  that  about  3.1  million  wetland  basins 
covering  about  28,490  km^  (Table  7)  were  in  the  prairie 
pothole  region  of  the  United  States.  Most  basins  were  in 
North  Dakota  (2.0  million),  and  more  were  in  South  Da¬ 
kota  (0.8  million)  than  in  Minnesota  (0.2  million).  Our 
sample  in  Montana  (48,800  wetland  basins)  represented 
only  a  small  area  of  that  state. 


Table  6.  Land  ownership  (%)  by  state  in  the  prairie  pothole  region  of  the  United  States  based  on  1986  data. 


Easement 


Minnesota 
Montana 
North  Dakota 
South  Dakota 
All  States 


^Includes  all  land  in  the  easement  tract.  Only  the  wetlands  are  under  easement. 

^Includes  government  land  not  owned  or  under  easement  by  the  U.S.  Fish  and  Wildlife  Service. 
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Table  7.  Estimated  number  of  basins,  area  of  wetland,  density  of  basins,  and  area  of  wetland/km  by  land-ownership  class  in  the  prairie  pothole  region  of  Minnesota, 

North  Dakota,  South  Dakota,  and  extreme  eastern  Montana,  based  on  photographs  from  the  early  1980s, 
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Hectares 


Fig.  5.  Cumulative  distribution  of  wet¬ 
land  basins  by  size  and  wetland-ba¬ 
sin  class  in  the  prairie  pothole  region 
of  the  United  States,  based  on  pho¬ 
tographs  from  the  early  1980s.  Wet¬ 
land-basin  classes  include  lake  ba¬ 
sins  (Lake),  semipermanent  wetland 
basins  (semipermanent),  seasonal 
wetland  basins  (seasonal),  and  tem¬ 
porary  wetland  basins  (temporary). 


2 

The  estimated  area  of  wetland  per  km  was  similar 
among  the  states  except  in  Montana  where  the  area  was 
about  half  of  that  in  other  states.  The  estimated  density  of 
wetland  basins  varied  among  states.  The  highest  density  was 
in  North  Dakota.  The  estimated  mean  wetland-basin  sizes 
were  2.7  ha  in  Minnesota,  1 .2  ha  in  Montana,  0.6  ha  in  North 
Dakota,  and  1.1  ha  in  South  Dakota.  The  distribution  of 
wetland  basins  by  size  class  was  highly  skewed  to  the 
smaller  sizes  (Fig.  5).  This  distribution  was  biased  because 
parts  of  some  wetland  basins  on  the  plot  boundary  included 
only  the  area  inside  the  plot.  The  bias  was  greatest  in  areas 
of  large  lakes.  The  estimated  wetland  area  per  km^  was 
greater  on  service  easements  than  on  private  land  and  was 


Easement  FWS  Private  All 


much  greater  on  service-owned  land  than  in  any  of  the  other 
landownership  classes  because  the  service  tended  to  buy 
large  wetland  basins,  especially  for  refuges.  Estimated  wet¬ 
land-basin  density  was  highest  on  service  easements  (Ta¬ 
ble  7)  because  easements  were  taken  in  areas  of  high  wet¬ 
land-basin  density  and  the  wetland  basins  were  small. 

The  distribution  of  wetland-basin  classes  varied  by  land- 
ownership  class  (Fig.  6).  The  distributions  of  wetland-basin 
classes  on  easements  and  on  private  land  were  similar, 
except  more  seasonal  wetland  basins  and  fewer  lakes  were 
on  easements  than  on  private  land.  Temporary  and  seasonal 
wetland  basins  were  greatly  underrepresented  and  lakes 
were  overrepresented  on  service-owned  land  most  of  which 

■  Lake 

B  Semipermanent 
B  Seasonal 
Q  Temporary 

Fig.  6.  Distribution  of  wetland-basin  classes  by 
landownership  class  in  the  prairie  pothole 
region  of  the  United  States,  based  on  photo¬ 
graphs  from  the  early  1980s.  Wetland-basin 
classes  include  lake  basins  (lake),  semiper¬ 
manent  wetland  basins  (semipermanent), 
seasonal  wetland  basins  (seasonal),  and  tem¬ 
porary  wetland  basins  (temporary).  Land- 
ownership  classes  include  land  owned  by  the 
U.S.  Fish  and  Wildlife  Service  (service)  and 
privately  owned  land  (private). 
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is  on  refuges,  which  were  historically  purchased  as  resting 
areas  for  migrating  birds  and  not  as  waterfowl  production 
areas. 

Upland  Habitat 

The  estimated  area  of  upland-habitat  classes  (Table  8) 
varied  by  state  and  landownership  class.  Cropland  was 
the  dominant  land-cover  class  on  private  land  and  on 
easements  in  all  four  states.  The  distribution  of  area  of 
land-cover  classes  on  easements  and  on  private  land  was 
similar.  The  small  amount  of  cropland  on  service-owned 
land  was  the  result  of  tillage  on  some  refuges  and  in 
waterfowl  production  areas  where  tillage  was  used  to 
rejuvenate  planted  cover  or  to  provide  food  plots.  These 
tilled  areas  appeared  as  cropland  on  aerial  photographs. 

Grassland  was  the  second  most  abundant  component 
of  the  landscape  in  all  states  and  its  abundance  increased 
from  east  to  west.  The  percent  of  grassland  was  lowest 
in  Minnesota  (6.7)  and  lower  in  North  Dakota  (16.3)  than 
in  South  Dakota  (26.0)  and  in  Montana  (33.9).  The 
distribution  of  grassland-wildlife  (Table  8)  was  partly 
an  artifact  of  our  definition.  No  grassland  was  shown  on 
service-owned  land  and  no  grassland-wildlife  was 
shown  on  easements  and  on  private  land.  The  occur¬ 
rences  of  hayland  and  other  land  were  minor  but  were 
important  duck  nesting  habitat  (Cowardin  et  al.  1985). 
Planted  cover  was  uncommon  except  on  service-owned 
land  where  it  was  a  major  habitat  component. 

Annual  Change  in  Ponds 

Pond  density  varied  from  4.4/km^  in  1987  to  0.8  in  1990 
(Fig.  7).  The  area  of  ponds  (ha  per  km^)  varied  from  7.2  in 
1987  to  2.7  in  1990  (Fig.  8).  Climate  in  the  prairie  pothole 
region  occurs  in  cycles  of  wetness  and  drought  (Kentaid  et  al 
1989).  Near-average  water  conditions  were  present  in  1987 
and  followed  by  drought  in  1988,  slight  recovery  in  1989, 
and  severe  drought  in  1 990.  The  highest  density  of  ponds 
was  on  easements,  especially  in  1987  (Fig.  7).  The  area  of 
ponds  per  km  was  greater  on  service-owned  land  than  in 
the  other  types  of  landownership  (Fig.  8).  Although  these 
large  basins  with  semipermanent  water  regimes  responded 
more  slowly  to  drought  conditions,  surface  water  declined 
steadily  throughout  the  period  (Fig.  8).  The  area  of  ponds 
per  km^  on  easements  and  on  private  land  fluctuated  simi¬ 
larly.  The  average  pond  size  (Fig.  9)  increased  slightly  as  the 
amount  of  surface  water  declined  because  small  wetland 
basins  dried  up  entirely.  The  average  pond  size  in  wetland 
basins  on  service-owned  lands  increased  during  1987-88, 
contracted  as  small  ponds  appeared  in  1989,  and  increased 
again  in  1990  as  the  small  ponds  again  disappeared.  The 
average  pond  size  was  less  in  1990  than  in  1988  because 
even  the  larger  ponds  were  becoming  smaller  because  of 
severe  drought. 


Duck  Populations  and  Production 

Breeding  Populations 

The  13  duck  species  reported  here  (Fig.  10)  represent 
39%  of  the  duck  species  that  breed  in  North  America  and 
73%  of  the  species  in  the  genera  Anas  and  Aythya  in  North 
America.  During  1987-90,  the  density  of  the  blue-winged 
teal  (3.4  pairs/km^)  was  higher  than  those  of  the  mallard, 
(2. 1 ),  gadwall  ( 1 .8),  northern  pintail  (0.8),  and  redhead  (0.8) 
(Fig.  1 1 ).  Density  was  consistently  highest  on  service-owned 
lands  primarily  because  of  the  large  areas  of  ponds. 

Annual  Change  in  Sizes  of  Breeding  Populations 

The  sizes  of  the  breeding  populations  of  the  five  most 
numerous  dabbling  ducks  declined  throughout  1 987-90  as 
drought  conditions  intensified  (Fig.  1 1).  The  declines  cor¬ 
responded  closely  to  loss  of  ponds  (Fig.  7)  and  pond  area 
(Fig.  8),  but  the  relation  of  pair  density  to  area  of  ponds 
differed  by  species  (Fig.  12).  We  expected  species  with  a 
high  degree  of  philopatry  and  possibly  weak  territoriality 
to  concentrate  on  ponds  as  the  number  of  ponds  decreased 
because  of  drought.  Slopes  of  the  linear  regressions  of  pair 
density  on  pond  density  (Fig.  13)  were:  -0.108  (gadwalls), 
-0.098  (mallards),  -0.071  (blue-winged  teals),  -0.025 
(northern  shovelers),  and  -0.0004  northern  pintails).  The 
ranking  is  similar  to  published  data  on  return  rates  (Ander¬ 
son  et  al.  1992).  Mallards  and  gadwalls  exhibit  strong 
philopatry  (Lakemoen  et  al  1990).  Our  ranking  of  blue¬ 
winged  teals  is  higher  and  of  northern  pintails  lower  than 
expected  according  to  the  literature.  The  comparison  with 
the  literature  can  only  be  approximate  because  of  consid¬ 
erable  spatial  and  temporal  variations  in  published  return 
rates,  variation  of  return  rates  by  age  and  sex,  and  rare 
correction  of  published  return  rates  of  mortality. 

Our  data  (Fig.  12)  may  suggest  that,  where  the  correla¬ 
tion  between  the  area  of  ponds  and  breeding  population  is 
low,  the  number  of  ducks  is  too  low  to  fill  the  available 
breeding  habitat.  However,  other  explanations  are  possi¬ 
ble,  and  we  found  no  data  that  support  a  depression  of 
northern  shoveler  populations.  Johnson  and  Shaffer  (1987) 
analyzed  data  from  annual  surveys  by  the  service  and 
concluded  that  estimated  mallard  population  sizes  no 
longer  parallel  estimated  pond  numbers.  Their  first  possi¬ 
ble  explanation  was  that  the  number  of  mallards  was  too 
low  to  fill  the  habitat. 

Change  in  the  total  number  of  pairs  per  km  by  year 
varied  among  landownership  classes  (Fig.  14).  The  high¬ 
est  density  was  on  service-owned  land  as  expected  be¬ 
cause  more  wetland  and  more  pond  area  are  on  these 
lands  than  on  land  in  the  other  landownership  classes. 
The  decrease  in  pairs  per  km^  during  1987-88  was  less 
severe  on  service-owned  lands  than  on  other  lands,  prob¬ 
ably  because  of  the  greater  amount  of  semipermanent 


14  Biological  Science  Report  2 


c 

o 

o 

CO 

O  CO 

0 

q 

0 

q 

0 

0 

q 

0 

q  q 

0 

On 

0 

ON 

1 

0 

0 

CO 

q 

0 

0 

to 

C3 

O 

2 

u 

o 

NO 

O 

d 

od 

d 

od 

d  NO 

d 

NO 

d 

in 

d  in 

d 

NO 

d 

NO 

> 

Cl 

d 

d 

NO 

d 

d 

NO 

^  NO 

r- 

r- 

r- 

CN 

CN 

NO 

NO 

a 

0 

0 

in 

in 

r- 

> 

p 

P 

d) 

E 

*0 

c 

ed 

o 

00 

q  q 

0 

0 

0 

0 

0  CN 

0 

q 

0 

q 

0  CN 

0 

0 

0 

0 

X 

0 

i 

95% 

0 

0 

q 

q 

0 

0 

§ 

d 

od 

d  od 

d 

d 

d 

d 

d  X 

d 

X 

d 

CN 

d  (N 

d 

d 

d 

d 

X 

d 

d 

in 

in 

d 

d 

u 

ON 

CNI 

CN 

NO 

so 

00 

00 

04 

CO 

CO 

>< 

CO 

CO 

d- 

d- 

*c 

q 

q 

(D 

oc 

X 

C 

03 

Cm 

O 

cd 

< 

fsl 

E 

o 

d 

NO 

0.0 

7.6 

0*0 

q 

od 

00 

q 

od 

0.0 

3.9 

0 

d 

q 

cd 

0.0 

q 

cd 

0.0 

3.6 

0 

d 

d; 

cd 

0 

d 

d 

cd 

ecj 

(N 

E 

0 

d 

0 

d 

q 

q 

0 

d 

0 

d 

Tt 

CO 

CO 

Tt 

ON 

On 

CN 

CN 

CN 

CN 

3 

'd- 

'd' 

NO 

NO 

< 

q 

q 

vJ 

NO 

o 

00  ^ 

q 

0 

q 

On 

0  0 

0 

q 

m 

0 

q  ^ 

q 

0 

O) 

q 

q 

q 

00 

0 

q 

C3 

d 

X  o 

06 

d 

CN 

ON 

ON  d 

NO 

X 

d 

d 

ON  cd 

d- 

d 

in 

d* 

U 

NO 

CN 

cd 

d 

d 

od 

Q 

u 

oo 

CO  00 

CO 

NO 

NO 

■d- 

0 

CO  d” 

NO 

in 

1— H 

> 

CO 

NO 

00 

m 

*0 

oo  00^ 

q 

q 

q 

CN 

q 

q 

q  0^ 

P 

q 

0^ 

0 

(N 

d 

c 

no"  no" 

r--" 

tX 

in" 

on" 

'd'" 

co" 

co"  rX 

00" 

cn" 

0 

ed 

CN 

CN  (N 

in 

SO 

'O 

in 

o 

c/5 

"c« 

(A 

ITD 

o 

(N  ON 

0 

q 

0  0 

Tj- 

q 

CO 

0 

CN  0 

0 

00 

00 

0 

c 

d; 

ON 

q 

r- 

0 

0 

«o 

00 

d 

d 

ON 

d 

cd 

ON 

d 

06 

od 

ON 

d 

d  NO 

CN 

d 

d* 

CN 

cd 

ON 

cd 

d 

in 

d 

d 

d 

c3 

OJQ 

ON 

CO 

ON  NO 

r- 

0 

0 

CO 

3 

in 

NO 

CN  00 

CO 

in 

'x 

in 

NO 

4» 

"Tf  to 

q 

q 

q 

q 

q 

d- 

q 

d; 

ON 

IM 

1M 

CN 

O 

O 

co"  CO 

X 

X 

co" 

■•st- 

00" 

cn" 

fX  0" 

so 

cn" 

ON 

0 

o 

cd 

d) 

M 

<{; 

CN 

E 

V£D 

o 

q  in 

q 

0 

q 

in 

q  q 

(N 

CN 

q 

0 

^  CN 

q  q 

d; 

0 

d- 

d; 

d 

q 

Cd 

E 

< 

CN 

c 

q 

CN 

q 

0 

q 

X 

X 

o 

CN 

od 

d 

iri 

cd  d 

tX 

d 

d- 

d 

ON  d- 

od 

d 

in 

cd 

c 

d 

NO 

ON 

in 

d 

cd 

NO 

NO  (N 

0 

CO 

Tt 

ON 

m 

in 

00 

r-  NO 

d- 

CO 

00 

X 

CN 

y-^ 

y-^ 

in 

Q 

<N 

q 

00^ 

q 

q 

q 

q 

q  q 

q 

q 

d^ 

CN 

y—f 

CO 

iz; 

in"  m" 

Tt" 

Tt" 

no" 

X 

cn" 

0"  co" 

tX 

rX 

m" 

CN 

CN  CN 

d- 

in 

a 

o 

C/D 

o 

cn 

q 

CN 

q 

q 

q 

q 

'd- 

NO 

q 

q 

q  q 

q 

CN 

q 

00 

q 

p 

d; 

q 

r-' 

CO 

C/D 

0) 

00 

X 

X 

CN  NO 

in 

in 

X 

in  od 

NO 

d 

d- 

rX 

NO  cd 

fX 

NO 

X 

ON 

U 

CN 

<N 

06 

iri 

d 

d 

ON 

r- 

00 

CO  CO 

'=3' 

CN 

r- 

r- 

d* 

00 

0  d- 

CO 

c^) 

H— 1 

2 

00 

vn 

CO 

in 

c 

U 

CO 

q 

q 

P  CO 

q 

■d-^ 

q 

d-^  p 

q 

d- 

q 

M 

»— 1 

q 

r-"  oo" 

X 

X 

CN 

tX 

ON 

in" 

ON  -d 

fX 

On" 

m" 

X 

in 

m" 

§ 

H! 

"O 

c 

in  in 

r- 

00 

d"  in 

00 

0 

CN 

Cd 

X 

Cm 

o3 

2 

u 

o 

d> 

p. 

NO 

(N 

H-H 

CN 

0 

q 

q 

q  q 

ON 

00 

d“ 

q 

00  d; 

q 

q 

d 

<3 

p 

d 

q 

q 

0 

0 

c 

d> 

o 

Cn)* 

d 

d 

in 

d 

CN 

NO 

Os 

On 

d 

cd 

od  od 

o^ 

•d 

X 

X 

M 

in 

cd 

CN 

vn 

od 

d 

d 

o 

«M 

o 

ir> 

o 

CO 

r-  in 

NO 

m 

VO  — 

0 

CO 

NO 

d  ^ 

NO 

CN 

so 

0 

On 

in 

CN 

CO 

CN 

Cm 

ON 

cn 

*'1 

q 

in  -- 

q 

ON 

q 

d;  q 

q 

CN 

q 

00 

in 

q 

*5) 

S2 

E 

in" 

in" 

0" 

On" 

0" 

co" 

X  m" 

d"" 

d-" 

0" 

0 

CO 

cd 

cd 

«n  in 

NO 

00 

d- 

r- 

ON 

Cd 

<3 

*0 

C/) 

< 

o 

q 

CN  OO 

Tf 

CN 

q 

q  c^^ 

q 

q 

NO 

r-  00 

q 

q 

q 

NO 

< 

CN 

q 

r- 

r-- 

0 

0 

Oh 

CN 

C 

d 

ui 

X  vo 

NO 

CN 

ON 

X 

tX 

cd 

in 

ri 

d 

rX  0 

cd 

in 

X 

d 

tX 

X 

NO 

CN 

d 

0 

Q 

3 

C 

'd- 

lO 

O  ON 

0 

'n 

ON  ^ 

ON 

0 

0 

in 

CN  00 

d- 

CN 

t-- 

d 

c 

X 

so 

CO 

CO 

d 

Oh 

M 

m 

q 

q  c^i 

q 

q 

q 

d^  q 

0 

CO 

p 

q 

q 

d^ 

d) 

;c 

04) 

OCD 

00" 

X 

CO 

in" 

d-" 

m"  On 

no" 

X 

od 

d" 

d-" 

O 

o 

m  in 

tM 

r- 

00 

d  d- 

00 

ON 

*c3 

X 

Oh 

Oh 

c 

00 

O) 

^  q 

00 

CN 

00 

q 

q 

q 

q 

NO 

q  q 

CN 

q 

00 

CO 

CO 

NO 

0 

0 

<D 

o 

to 

X 

CN 

d 

d 

'd- 

in 

od  CN 

X 

NO 

cd 

d 

cd  in 

CN 

cd 

On 

d 

in 

X 

d 

ON 

d 

d 

'O 

r  ) 

NO 

NO 

CO 

On 

On 

VO 

CN 

CN 

CN 

u  1 

CN 

On 

0 

r-; 

q 

q 

00 

CN 

q 

q 

q^ 

.s 

d) 

C/D 

C 

CO  CO 

'd'" 

-d-" 

X  X 

ON 

ON 

(N 

CN 

C/D 

C/D 

03 

<3  , 

X  i 

c 

CO 

^  q 

CO 

0 

00 

, 

0  CN 

0 

q 

q 

q  q 

CN 

q 

ON 

q 

1 

q 

0 

q 

q 

0 

0 

13 

ci 

ri 

d 

d  in 

d 

d 

CN 

cd 

NO  d 

NO 

X 

od 

d 

od  On 

X 

ON 

X 

cd 

in 

d 

vn 

X 

0 

d 

C 

00  oo 

NO 

NO 

On 

ON 

0 

(N 

CO  0 

r- 

NO 

X 

ON 

<N 

d" 

_Oh 

iS 

X 

On 

q^  «p 

00 

0^ 

d;  q 

m 

q 

q 

q 

c 

o 

X  X 

X 

y—i  mH 

d-" 

d-" 

0 

mh 

»— 1 

C/D 

Ui 

<u 

o 

cd 

d) 

u. 

CN 

cn 

1 

00 

NO 

. 

OO 

q  q 

q 

q 

q 

CO 

2  ^ 

q 

q 

p 

NO 

cd 

E 

< 

CN 

c 

q 

q 

q 

q 

0 

0 

O 

< 

£ 

X 

vd  X 

d 

d 

cd 

X  X 

cd 

X 

X 

d 

X  (N 

cd 

CN 

d- 

d 

in 

CN 

cX 

in 

CD 

d 

ON  o 

NO 

NO 

0  NO 

00 

in 

d 

c 

m 

r- 

o 

q 

*M 

1— H 

h~h 

q 

q 

q  q^ 

1— 1 

q 

0^ 

X 

q 

0^ 

-o 

o 

c3 

CN  CN 

cn" 

CN 

so 

tX 

CN 

cn" 

X) 

C/D 

d) 

d> 

d> 

0 

0 

CL) 

<3 

s 

o 

0 

u 

0 

u 

0 

.M 

*E 

> 

> 

'E 

'> 

’t 

X 

d> 

oS 

0 

<u 

3 

c:3 

{/3 

c/3 

t/3 

00 

c/3 

C/3 

C/3 

X 

di 

d) 

d) 

0) 

«D 

.0 

,<3 

X 

nm 

X 

X 

iM 

nm 

nm 

Oh 

X 

'X 

X 

X 

X 

X 

O 

X 

2 

2 

2 

2 

2 

2 

2 

2 

Oh 

3 

0^ 

1 

1 

1 

3 

i 

Cm 

C 

73 

13 

73 

73 

X) 

c 

X3 

T3 

Q 

O 

C 

C 

c 

c 

a 

C 

03 

d) 

T3 

o 

cd 

cd 

cd 

Cd 

cd 

0 

cd 

cd 

C 

d) 

E 

d) 

c/5 

X 

C/D 

00 

cd  33 

>  i2 

s  .’i 

E 

2^  to 

cd 

> 

2 

S  '§ 

£  X 

d> 

iA  c/3 

2 

cd 

> 

2 

c 

<L> 

E 

0 

\A 

X 

<j0 

X 

od 

M 

2 

>  B 

c 

<D 

E 

(D 

lA 

X 

on 

cd 

M 

2 

> 

2 

S  "§ 

£  X 
on  cd 

M 

2 

> 

*2 

c 

<3 

E 

<3 

tA 

X 

<Z5 

X 

C/3 

d> 

cd 

’C  0 

cd 

'C 

*C 

0 

Cd 

'C  0 

C3 

'C 

0 

cd 

0 

Cd 

P 

W 

D 

W 

D 

X 

H 

X  D 

X 

X 

X 

X  H 

X 

D 

X 

H 

X 

D 

X 

H 

X 

X 

M 

C/) 

cd 

w 

M 

0 

0 

00 

a> 

X 

CQ 

H 

State 

Minnesota 

Montana 

0 

’3 

tn 

0 

2 

P 

X 

3 

0 

00 

All  States 

State 

Minnesota 

Montana 

Private  1.9  0.0  5,7  62.6  25.3  100.0  3.1  0.0  8.6  103.7  34.1  173.2 

Total  1.9  0.0  5.7  63.0  25.7  100.2  6.9  0.0  19.9  103.7  34.1  173.2 


North  Dakota  Easement  582.3  382.1  782.5  398.5  339.2  457.7  103.3  0.0  218.9  14.0  2.5  25.5 

U.S.  Fish  and  Wildlife  Service  1.5  0.0  4.7  113.2  0.0  231.4  202.7  105.7  299.7  2.6  0.4  4.8 

Private  2,490.2  1,850.4  3,130.2  3,849.2  2,122.1  5,576.3  293.5  157.0  430.0  1,881.0  1,625.3  2,136.8 

Total  3,074.1  2,365.2  3,782.9  4,360.8  2,492.9  6,228.7  599.4  369.2  829.7  1,897.6  1,636.4  2,158.8 


Evaluations  of  Duck  Habitat  and  Estimation  of  Duck  Population  Sizes  1 5 


CO 

o 

r- 

CM 

CO 

oo 

q 

q 

c4 

d 

d 

d 

wo 

WO 

CM 

CO 

NO 

q 

q 

q 

wo" 

wo" 

o 

o 

q 

1— 

Tj- 

CM 

q 

d 

d 

wo 

CO 

d 

CM 

00 

8 

CO 

CO 

q 

r- 

o 

o 

oo 

q 

NO 

oo 

q 

d 

d 

o\ 

as 

CM 

CM 

CNJ 

CM 

CO 

WO 

CO 

q 

q 

q 

q 

wo 

— 

CM 

r- 

q 

w^ 

wo 

d 

d 

CO 

CO 

d 

Tt 

d 

On 

(N 

ON 

CM 

r-- 

so 

CM 

NO 

On 

q 

O 

o 

rf 

CO 

q 

q 

d 

d 

d 

d 

d 

wo 

d 

CO 

CN 

wo 

CO 

r- 

Tj- 

CM 

CO 

CO 

ON 

CO 

o 

r- 

oo 

CM 

q 

00 

CO 

d 

Tf 

d 

wo 

d 

r- 

CN 

CO 

NO 

On 

CO 

C3N 

CN 

CO 

NO 

q 

WO 

CO 

oo 

q 

q 

d 

oi 

d 

d 

CO 

od 

00 

1-^ 

Tt 

ON 

00 

CO 

1 

q 

NO 

NO 

CM 

CM 

<— 1 

cs 

cm" 

cm" 

co" 

(S 

00 

CN 

q 

q 

CN 

S 

00 

On 

1—^ 

o^ 

CM 

CO 

ON 

CO 

q 

WO 

q 

q 

oo" 

oo 

CO 

o 

CN 

ON 

r- 

q 

d 

CN 

ON 

WO 

q6 

d 

wo 

CO 

CM 

NO 

o 

w^ 

r- 

CO 

O 

so 

CM 

q 

CM 

CM 

o" 

q 

q 

f" 

CM 

r- 

1—1 

00 

CM 

'=t 

d 

d 

wo 

00 

00 

CO 

t-- 

CM 

r- 

q^ 

(M 

Tt 

wo" 

t-iT 

ON 

d 

r- 1 

CO 

o 

q 

CO 

Tt 

wo 

d 

wS 

d 

8 

ON 

On 

wo 

WO 

00 

CM^ 

NO 

q 

cm" 

d 

r-" 

wo 

O) 

q 

ON 

q 

q 

Td- 

wo 

00 

CO 

00 

d 

wo 

oo 

00 

CM 

00 

s 

CO 

r- 

ON 

q 

CO 

co" 

r-" 

oo" 

0) 

CL> 

u 

2 

■r 

’> 

C 

(L> 

u 

o 

00 

CO 

1<C) 

2 

1 

•o 

C 

c 

cd 

cd 

c 

£ 

t/3 

S 

E  P- 

a 

Cd 

o 

(/3 

00 

> 

Id 

(U 

(XI 

CO 

> 

"rt 

cd 

‘C 

o 

cd 

’C 

O 

u 

D 

o- 

H 

W 

H 

cd 

o 

c/3 

Q 

Oj 

Cd 

•S 

s 

o 

00 

< 

Area  of  ponds/km^  (ha) 


1 6  Biological  Science  Report  2 


1987  1988  1989  1990 

Year 


Fig.  7.  Density  of  ponds  by  land- 
ownership  class  during 
1987-90  in  the  prairie  pot¬ 
hole  region  of  the  United 
States,  based  on  aerial 
videography.  Landownership 
classes  include  easements 
(easement)  where  the  U.S. 
Fish  and  Wildlife  Service  has 
easements  on  wetlands  to 
prevent  draining,  filling,  or 
leveling;  land  owned  by  the 
U.S.  Fish  and  Wildlife  Serv¬ 
ice  (service);  privately 
owned  land  (Private);  and  all 
landownerships  combined 
(all).  Error  bars  =  1  standard 
error. 
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Fig.  8.  Area  of  ponds  per  km^ 
by  landownership  class, 
during  1 987-90  in  the  prai¬ 
rie  pothole  region  of  the 
United  States,  based  on 
aerial  videography.  Land- 
landownership  classes  in¬ 
clude  easements  (ease¬ 
ment)  where  the  U.S.  Fish 
and  Wildlife  Service  has 
easements  on  wetlands  to 
prevent  draining,  filling,  or 
leveling;  land  owned  by  the 
U.S.  Fish  and  Wildlife 
Service  (service);  privately 
owned  land  (private);  and 
all  landownerships  com¬ 
bined  (All).  Error  bars  -  1 
standard  error. 
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Fig.  9.  Average  size  of  ponds 
by  landownership  class 
during  1987-90  in  the  prai¬ 
rie  pothole  region  of 
the  United  States,  based  on 
aerial  videography.  Land- 
landownership  classes  in¬ 
clude  easements  (ease¬ 
ment)  where  the  U.S.  Fish 
and  Wildlife  Service  has 
easements  on  wetlands  to 
prevent  draining,  filling,  or 
leveling;  land  owned  by  the 
U.S.  Fish  and  Wildlife 
Service  (service);  privately 
owned  land  (private);  and 
all  landownerships  com¬ 
bined  (all).  1987-90. 
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Fig.  10.  Average  density  of  breeding  ducks  (Anatinae)  by  landownership  class  during  1987-90  in  the  prairie  pothole  region  of  the 
United  States,  based  on  model  projections  from  aerial  videography  and  ground  counts.  Duck  species  include  mallards  {Anas 
platyrhynchos),  gadwalls  {Anas  strepera)^  American  wigeons  {Anas  americana),  green-winged  teals  {Anas  crecca),  blue-winged 
teals  {Anas  cJiscors),  northern  shovelers  {Anas  clypeata),  northern  pintails  {Anas  acuta),  wood  ducks  {Aix  sponsa),  redheads 
{Aythya  americana),  canvasbacks  {Aythya  valisnneria),  lesser  scaups  {Aythya  affinis),  ring-necked  ducks  {Aythya  coUaris),  and 
niddy  ducks  {Oxyura  jamaicensis). 
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Fig.  11.  Relation  of  density  of 
blue-winged  teal  {Anas  clis- 
cors),  gad  wall  {Anas  strepera), 
mallard  {Anas  pkityrhynchos), 
northern  shoveler  {Anas 
clypeata),  and  northern  pintail 
{Anas  acuta)  pairs  to  pond  size 
during  1987-90  in  the  prairie 
pothole  region  of  the  United 
States.  Error  bars  =  1  standard 


Year 


Fig.  12.  Relation  between  density  of  blue- winged 
teal  {Anas  discors),  gadwall  {Anas  streperd),  mal¬ 
lard  {Anas  platyrhynchos),  northern  shoveler 
{Anas  clypeata),  and  northern  pintail  {Anas  acuta) 
pairs  to  area  of  ponds  during  1987-90  in  the 
prairie  pothole  region  of  the  United  States. 
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Fig.  13.  Relation  between  average  number  of  pairs  of 
blue-winged  teals  {Anas  discors),  gadwalls  {Anas 
strepera),  mallards  {Anas  platyrhynchos),  north¬ 
ern  shovelers  {Anas  clypeata),  and  northern  pin¬ 
tails  {Anas  acuta)  per  pond  and  pond  density 
during  1987-90  in  the  prairie  pothole  region  of 
the  United  States. 


Fig.  14.  Relation  between  the  av¬ 
erage  number  of  duck  (Anati- 
nae)  pairs  per  km^  and  year  by 
landownership  class  during 
1987-90  in  the  prairie  pothole 
region  of  the  United  States. 
Landownership  classes  in¬ 
clude  easements  (easement) 
where  the  U.S.  Fish  and  Wild¬ 
life  Service  has  easements  on 
wetlands  to  prevent  draining, 
filling,  or  leveling;  land  owned 
by  the  U.S.  Fish  and  Wildlife 
Service  (service);  privately 
owned  land  (private). 
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wetland  basins  on  service-owned  land.  Birds  that  return 
to  a  landscape  in  a  drought,  where  the  less  permanent 
ponds  were  dry,  probably  concentrated  on  large  semiper¬ 
manent  wetland  basins  like  those  on  service-owned  land. 
The  decline  of  the  number  of  pairs  was  more  apparent  on 
easements  than  on  private  lands.  More  wetland  basins 
are  on  easements  than  on  private  lands,  but  not  as  many 
large  semipermanent  and  permanent  wetland  basins  are 
on  easements  as  on  private  lands  (Fig.  8). 

Duck  Production 

The  number  of  recruits  is  the  product  of  size  of  breed¬ 
ing  population  and  the  recruitment  rate  (Cowardin  and 
Johnson  1979).  Drought  has  a  negative  effect  on  both 
(Cowardin  et  al.  1985,  Johnson  and  Grier  1988).  We 
point  out  that  our  estimated  recruitment  rate  was  more 
dependent  on  model  prediction  than  on  observation  and 
was  highly  influenced  by  the  underlying  assumptions  of 
the  model.  The  estimated  density  of  the  recruited  duck¬ 
lings  (Fig.  15)  followed  the  same  general  pattern  as  the 
sizes  of  the  breeding  populations  among  the  five  species 
(Figs.  10  and  15)  and  among  years  (Figs.  1 1  and  16)  of 
the  five  species  for  which  production  was  estimated.  The 
estimated  recruitment  rates  varied  among  species  and 
among  years  (Fig.  17).  The  rates  were  highest  in  blue¬ 
winged  teals  and  gadwalls  and  lowest  in  mallards  and 
northern  pintails.  The  annual  variation  in  recruitment 
rates  (Fig.  17)  resulted  from  variation  in  a  (a  measure  of 
nesting  intensity;  Table  9).  Our  estimates  of  a  had  a 
major  effect  on  our  estimates  of  hen  success  (Equa¬ 


tion  5).  The  estimated  clutch  success  in  stable  popula¬ 
tions  (Cowardin  et  al.  1 985  and  Klett  et  al.  1988)  is  lower 
than  those  presented  in  Table  10.  However,  the  estimates 
presented  in  those  papers  were  based  on  the  assumption 
that  a  equaled  1.  For  the  low  a  values  in  this  study, 
higher  clutch  success  is  required  for  recruitment  rates  of 
a  stable  population.  For  mallards,  a  hen  success  of  3 1  %, 
a  summer  survival  of  0.74,  and  an  average  brood  size  of 
4.9  (Cowardin  et  al.  1985)  equate  to  a  recruitment  rate 
of  0.56  in  a  stable  population.  This  is  well  above  our 
estimated  recruitment  rate  in  mallards  (Fig.  17). 

We  did  not  have  sufficient  data  for  estimating  clutch 
success  by  year.  Greenwood  et  al.  (1995)  showed  that 
clutch  success  in  the  prairies  of  Canada  is  depressed  by 
drought.  If  we  had  estimates  of  clutch  success  by  year, 
the  variation  in  our  estimated  recruitment  rate  by  year 
(Fig.  17)  would  probably  have  been  greater. 

Our  model  predicted  that  the  densities  of  mallard, 
gadwall,  blue-winged  teal,  and  northern  shoveler  re¬ 
cruits  were  greater  on  service-owned  land  than  on  ease¬ 
ments  and  private  land  (Fig.  16).  Northern  pintail  re¬ 
cruits  were  produced  at  almost  equal  rates  in  the 
landownership  classes.  The  number  of  successful  nest¬ 
ing  attempts  per  km^  varied  by  landownership  class 
(Fig.  18)  because  of  the  differences  in  clutch  success 
(Table  10)  and  the  relative  preferences  of  the  different 
species  for  the  various  nesting  habitats  (Klett  et  al.  1988, 
Table  2)  on  those  lands.  Klett  et  al.  (1988)  found  a  higher 
preference  for  cropland  in  northern  pintails  than  in 
other  duck  species.  The  densities  of  successful  nesting 


Fig.  15.  Number  of  produced  recruits 
per  km^  by  mallards  {Anas 
platyrhynchos) ,  gadwalls  {Anas 
strepera),  blue-winged  teals 
{Anas  discors),  northern  shovel- 
ers  {Anas  clypeata),  and  north¬ 
ern  pintails  {Anas  acuta)  during 
1987-90  in  the  prairie  pothole  re¬ 
gion  of  the  United  States,  based  on 
model  projections  from  aerial 
videography  and  ground  counts. 
Landownership  classes  include 
easements  (easement)  where  the 
U.S.  Fish  and  Wildlife  Service  has 
easements  on  wetlands  to  prevent 
draining,  filling,  or  leveling;  land 
owned  by  the  U.S.  Fish  and  Wild¬ 
life  Service  (service);  privately 
owned  land  (private);  and  all  land- 
ownerships  combined  (all). 
1987-90.  Error  bars  -  1  standard 


error. 
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Fig.  16.  Number  of  produced  re- 
cruits  per  mi^  of  blue¬ 
winged  teals  (Anas  discors), 
gadwalls  (Anas strepera),m2i\- 
lards  (Anas  platyrhynchos)^ 
northern  shovelers  (Anas 
clypeata),  and  northern  pin¬ 
tails  (Anas  acuta)  during 
1 987-90  in  the  prairie  pothole 
region  of  the  United  States, 
based  on  model  projections 
from  aerial  videography  and 
ground  counts.  Error  bars  =  1 
standard  error. 
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Fig.  17.  Recruitment  rate  of  five 
species  of  dabbling  ducks  dur¬ 
ing  1 987-90  in  the  prairie  pot¬ 
hole  region  of  the  United  States, 
based  on  model  projections  from 
aerial  videography  and  ground 
counts.  Species  include  mallards 
(Anas  platyrhynchos),  gadwalls 
(Anas  strepera),  blue-winged 
teals  (Anas  discors),  northern 
shovelers  (Anas  clypeata),  and 
northern  pintails  (Anas  acuta) 
and  all  species  combined. 


Year 
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Table  9.  Estimates  of  a,  a  measure  of  nesting  intensity,  by  wetland  management  district  and  year  during  1987-90  in 

the  prairie  pothole  region  of  the  United  States.^ 


District 

1987 

1988 

1989 

1990 

Arrowwood 

0.962 

0.776 

0.711 

0.668 

Audubon 

0.858 

0.721 

0.655 

0.640 

Crosby-Lostwood 

0.736 

0.655 

0.729 

0.629 

Devils  Lake 

0.745 

0.650 

0.656 

0.642 

J.  Clark  Salyer 

0.744 

0.658 

0.654 

0.630 

Kiilm 

0.879 

0.673 

0.759 

0.657 

Long  Lake 

0.920 

0.701 

0.678 

0.640 

Medicine  Lake 

0.820 

0.669 

0.844 

0.697 

Tewaukon 

0.741 

0.751 

0.766 

0.676 

Waubay 

0.803 

0.684 

0.844 

0.745 

is  a  linear  function  of  the  percentage  of  wetland  basins  containing  water  in  each  year. 


Table  10.  Estimates  of  average  clutch  success  by  species  and  landownership  class  during  1 987-90  in  the  prairie  pothole 

region  of  the  United  States. 


Species 

Ownership 

Easement 

FWS 

Private 

All  Owners 

Mallard 

17 

22 

17 

17 

Gadwall 

22 

27 

24 

24 

Blue-winged  teal 

22 

23 

23 

23 

Northern  shoveler 

21 

21 

21 

21 

Northern  pintail 

20 

21 

21 

21 

All  species 

20 

24 

21 

21 
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Fig.  18,  Average  number  of  suc¬ 
cessful  nesting  attempts  per 
km^  by  mallards  (Anas 
platyrhynchos),  gadwalls  (Anas 
strepera),  blue- winged  teals 
(Anas  discors),  northern  shov¬ 
el  ers  (Anas  clypeata),  and 
northern  pintails  (Anas  acuta) 
by  landownership  class  in  the 
prairie  pothole  region  of  the 
United  States,  based  on  model 
projections  from  aerial  videogra- 
phy  and  ground  counts.  Land- 
ownership  classes  include 
easements  (easement)  where 
the  U.S.  Fish  and  Wildlife 
Service  has  easements  on  wet¬ 
lands  to  prevent  draining,  fill¬ 
ing,  or  leveling;  land  owned  by 
the  U.S.  Fish  and  Wildlife 
Service  (service);  and  privately 
owned  land  (private). 
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attempts  per  km^  in  mallards,  gadwalls,  and  blue-winged 
teals  were  higher  on  service-owned  lands  than  on  other 
lands.  The  number  of  successful  nesting  attempts  per 
km^  by  northern  pintails  and  northern  shovelers  was 
similar  among  the  landownershipdasses. 

Conclusions  and  Management 
Implications 

Wetland  and  Duck  Management 

The  Effect  of  Service-owned  Land 

Sidle  (1983)  described  how  alarm  over  declining  duck 
populations  and  destruction  of  prairie-pothole  wetland  habi¬ 
tat  led  to  the  service’s  small  wetland  acquisition  program 
(DeBates  1967).  Our  data  demonstrated  how  waterfowl 
production  areas  and  service  easements  acquired  under  this 
program  and  national  wildlife  refuges  are  important.  These 
lands  not  only  prevent  further  destruction  of  wetland  but  also 
provide  areas  where  ducks  can  continue  to  reproduce.  If 
these  lands  can  be  managed  so  that  recruitment  rate  more 
than  compensates  for  annual  mortality,  the  managed  lands 
may  slow  or  reverse  the  declines  of  some  duck  populations. 
In  10  waterfowl  management  districts,  service-owned  land 
was  only  1 .3%  of  the  land  surface,  but  2.5%  of  the  five  most 
common  dabbling  ducks  were  produced  on  it.  On  easements, 
including  the  wetlands  under  easement  and  the  private  land 
in  the  easement  tract,  19.6%  of  the  ducks  were  produced  on 
14.1%  of  the  land  surface.  Private  land  on  the  other  hand 
produced  77.9%  of  the  ducks  on  84.6%  of  the  land  surface. 
These  estimates  do  not  include  the  contribution  of  service- 
owned  land  and  easements  where  the  preserved  wetlands 
contribute  to  duck  production  from  surrounding  private 
uplands.  Thus,  lands  owned  in  fee  and  easements  taken  on 
wetlands  increase  duck  production;  however,  most  ducks 
are  produced  on  private  lands. 

The  people  who  initiated  the  small  wetland-acquisition 
realized  that  bigger  and  wetter  wetland  is  not  necessarily 
better  for  ducks.  Kaminski  and  Weller  (1992)  summarized 
results  from  several  studies  that  revealed  the  importance  of 
seasonal  wetlands  to  dabbling  ducks.  These  wetlands  are 
generally  small.  The  cumulative  frequency  distribution  of 
wetland  basins  by  size  class  in  the  prairie  pothole  region 
(Fig.  5)  has  important  management  and  conservation  impli¬ 
cations.  Restricting  regulatory  responsibility  or  protection  to 
larger  basins  would  leave  most  of  the  prairie  potholes  un¬ 
protected.  For  example,  78%  of  the  basins  are  less  than  0.41 
ha  and  most  have  temporary  or  seasonal  water  regimes. 

Wetland  easements  are  the  main  source  of  protection  for 
wetland  basins  with  temporary  and  seasonal  water  regimes 
because  these  classes  are  poorly  represented  on  lands  owned 
in  fee  title  (Fig.  6).  The  large  lakes  on  national  wildlife 


refuges  may  serve  as  resting  areas  for  migrating  waterfowl, 
but  their  contribution  to  duck  production  is  minor. 

Drought  Effects 

Our  estimates  from  the  remote-sensing-based  system 
during  the  first  4  years  clearly  supported  the  well  known 
fact  that  drought  depresses  duck  populations  and  produc¬ 
tion  in  the  prairie.  Data  gathered  during  this  period  re¬ 
vealed  only  the  decline  in  breeding  populations  as  drought 
conditions  increased  during  1987-90.  Data  gathered  in 
future  years  may  document  the  increase  in  populations  and 
production  as  water  conditions  improve  on  the  prairies. 
Although  the  large  wetland  basins  on  service-owned  lands 
contribute  little  to  duck  production,  our  data  suggested 
that,  because  they  may  remain  wet  during  drought,  large 
wetland  basins  with  permanent  water  regimes  may  modify 
drought  effects. 

Species  Effects 

Duck  management  is  often  based  on  mallards  because 
data  are  more  comprehensive  on  mallards  than  on  other 
duck  species.  The  assumption  is  that  management  that  is 
good  for  mallards  is  good  for  other  ducks.  The  assumption 
is  probably  valid  for  some  dabbling  ducks,  but  managers 
should  be  aware  of  species  differences.  Our  analysis 
agreed  with  published  information  on  philopatry 
(Lokemoen  et  al.  1990),  which  suggests  that  homing  to 
natal  areas  is  stronger  in  mallards  and  gadwalls  than  in  the 
other  species  we  studied.  This  means  that  intensive  man¬ 
agement  to  increase  local  recruitment  would  be  more 
beneficial  to  these  species  than  to  those  that  do  not  home 
strongly.  ^ 

The  analysis  of  successful  nesting  attempts  per  km“  by 
landownership  class  (Fig.  18)  revealed  additional  important 
species  differences  from  the  interaction  of  preferences  for 
nesting  covers,  clutch  success  in  the  various  covers,  and  the 
types  of  cover  in  the  landownership  classes.  Accordingly, 
management  of  uplands  to  benefit  mallards  and  gadwalls  on 
service-owned  lands  may  not  be  beneficial  to  northern  pin¬ 
tails  and  northern  shovelers.  Upland  management  probably 
has  limited  value  for  diving  ducks,  most  of  which  nest  in 
wetlands.  Management  of  service-owned  lands  and  ease¬ 
ments  to  preserve  wetland  may  have  a  beneficial  effect  by 
providing  nesting  areas  for  diving  ducks  and  for  some 
dabbling  ducks,  especially  in  dry  years. 

System  Evaluation 

Possible  Biases 

The  estimates  derived  from  the  remote-sensing-based 
system  may  have  several  biases.  Testing  for  these  biases 
requires  special  studies  that  were  beyond  the  scope  of  our 
evaluation.  Typical  difficulties  in  estimating  duck  numbers 
and  production  over  vast  areas  may  preclude  objective  veri- 
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fication.  Our  purpose  here  is  to  alert  the  reader  to  some  of 
the  more  important  possible  sources  of  errors. 

The  remote-sensing-based  systehi  is  based  on  data  from 
remote  sensing  for  the  identification  of  wetlands  and  up¬ 
lands  from  aerial  photographs  and  on  data  from  aerial 
video  for  measuring  annual  change  in  the  number  of  ponds 
and  area  of  water.  We  made  the  assumption  that  mapping 
by  the  National  Wetland  Inventory  was  without  errors  of 
omission  or  commission.  Although  not  tested,  our  experi¬ 
ence  in  this  and  numerous  other  studies  in  the  prairie 
pothole  region  revealed  that  such  errors  are  few.  Errors 
that  we  detected  were  usually  errors  of  omission  of  small 
temporary  wetland  basins.  These  errors  lead  to  underesti¬ 
mation  of  breeding-population  sizes.  Data  on  the  number 
of  ponds  and  on  the  area  of  water  from  video  (sometimes 
poor  quality)  may  contain  errors  from  the  omission  of 
numbers  of  ponds  (L.  L.  Strong,  Northern  Prairie  Science 
Center,  National  Biological  Service,  Jamestown,  N.  Dak., 
unpublished  data).  These  errors  cause  the  underestimation 
of  breeding-population  sizes.  They  also  cause  the  under¬ 
estimation  of  a  and,  therefore,  recruitment  rate.  We  also 
point  out  that  our  method  for  estimating  a  relies  on  data 
from  only  4  years  in  one  study.  We  also  make  the  assump¬ 
tion  that  the  form  of  the  regression  curves  for  estimating 
the  number  of  breeding  pairs  is  the  same  among  years  and 
areas.  This  assumption  was  not  examined.  Adjustment  of 
regressions  by  y  may  have  introduced  errors  that  we  are 
currently  unable  to  evaluate. 

For  the  estimates  of  duck  production,  we  used  Conser¬ 
vation  Reserve  Program  cover  estimates  from  1990  and 
applied  them  to  all  years.  Conservation  Reserve  Program 
contracts  were  taken  during  1987-90.  Conservation  Re¬ 
serve  Program  cover  is  attractive  to  ducks  with  relatively 
high  clutch  success  (Kantrud  1993,  Reynolds  et  al.  1994); 
therefore,  our  estimates  of  production  in  years  prior  to 
1990  may  be  high. 

Ail  systems  for  the  estimation  of  breeding-population 
sizes  of  ducks  are  subject  to  errors  of  biological  interpre¬ 
tation  such  as  determining  whether  an  observed  pair  of 
ducks  represents  a  resident  or  migrant  pair  and  whether 
observed  social  groups  such  as  lone  males  represent  breed¬ 
ing  pairs  (Cowardin  and  Blohm  1992).  We  made  no  at¬ 
tempt  to  solve  these  problems  but  attempted  to  use  meth¬ 
ods  that  correspond  to  those  used  in  other  surveys  so  that 
the  estimates  are  comparable. 

Sampling  Errors 

Although  our  overall  sample  size  (Table  1 )  was  large, 
samples  in  some  wetland-management-district-land- 
ownership  strata  were  small,  and  variation  in  habitats 
among  plots  was  great.  Confidence  limits  on  most  esti¬ 
mates  of  habitat  parameters  were  large.  Cost  probably 
prohibits  large  increases  of  sample  size,  although  addi¬ 


tional  plots  are  currently  added  to  the  sample  where  sam¬ 
ple  sizes  are  minimal.  At  present,  wetland  mapping  and 
digitizing  in  the  entire  prairie  pothole  region  are  progress¬ 
ing  rapidly,  and  technology  is  available  for  adding  uplands 
to  the  data.  When  data  from  the  entire  prairie  pothole 
region  become  available,  sampling  will  no  longer  be  nec¬ 
essary  for  estimating  numbers  of  wetland  basins  and  wet¬ 
land  area.  Our  confidence  limits  for  duck-population  sizes 
and  duck-production  parameters  reflect  only  the  variation 
associated  with  our  sample  of  10.4-km^  plots.  For  exam¬ 
ple,  the  confidence  limits  do  not  reflect  the  variance  in  the 
number  of  ducks  or  recruitment  in  individual  plots.  Thus, 
the  confidence  limits  for  these  parameters  are  actually 
narrower  than  they  should  be. 

Suggested  Future  Modifications 

The  remote-sensing-based  system,  although  opera¬ 
tional,  should  be  considered  a  prototype  that  can  be  im¬ 
proved  by  various  modifications  and  by  the  addition  of 
new  data.  One  simple  modification  would  be  a  complete 
inventory  of  wetlands  as  described  in  the  previous  section. 
Annual  estimation  of  the  number  of  ponds  will  probably 
have  to  remain  sample  based  because  of  cost  and  logistic 
problems  in  obtaining  complete  remote-sensing  data  over 
such  a  large  area  each  spring.  Satellite  data  with  sufficient 
resolution  to  delineate  the  small  ponds  that  ducks  use  may 
become  available  at  a  reasonable  cost  in  the  future,  but 
obtaining  data  in  the  required  narrow  time  frame  will 
remain  a  problem.  Advances  in  aerial  videography  tech¬ 
niques  have  been  made  since  the  remote-sensing-based 
system  was  initiated.  We  recommend  that  new  methods  be 
evaluated  to  increase  the  accuracy  of  the  annual  estimates 
of  the  number  of  ponds.  We  used  ratio-estimation  methods 
involving  the  amount  of  land  area  in  each  landownership 
class  to  improve  our  estimates  of  several  parameters.  Ratio 
estimates  involving  area  of  wetland  or  numbers  of  basins 
in  each  landownership  class  could  further  improve  esti¬ 
mates  of  parameters  such  as  area  and  number  of  ponds  or 
numbers  of  breeding  pairs.  Such  improvements  will  be 
possible  when  the  inventory  of  wetlands  is  complete. 

As  use  of  the  remote-sensing-based  system  continues, 
new  data  become  available  for  improving  the  regression 
equations  for  estimating  numbers  of  breeding  pairs  from 
pond  areas.  Furthermore,  the  classification  of  wetland 
basins  should  be  modified  to  better  represent  major  differ¬ 
ences  in  duck  use.  For  example,  the  current  class  lakes 
includes  many  types  of  lakes  that  receive  different  use  by 
ducks.  The  data  should  be  examined  to  determine  whether 
enough  information  exists  to  construct  regressions  for 
various  kinds  of  lakes.  In  addition,  the  remote-sensing- 
based  system  was  designed  primarily  for  estimating  num¬ 
bers  of  pairs  on  discrete  basins.  The  remote-sensing-based 
system  does  not  work  well  for  estimating  numbers  of  pairs 
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on  large  lakes  and  in  riverine  habitat.  Studies  should  be 
conducted  to  determine  whether  shoreline  segments  of 
lakes  and  reaches  of  rivers  are  preferable  sampling  units 
in  these  habitats.  Such  a  change  requires  gathering  data 
from  these  units  to  develop  appropriate  regression  models 
for  the  remote-sensing-based  system. 

When  we  designed  the  remote-sensing-based  system, 
we  could  not  link  vector  data  from  the  National  Wetland 
Inventory  maps  to  the  raster  data  from  aerial  videography. 
Therefore,  ground  crews  estimated  the  proportion  of  a 
wetland  basin  containing  water.  These  estimates  were  used 
to  calculate  y.  The  link  between  the  two  data  sets  is  now 
available,  and  we  recommend  that  aerial  videography  (or 
possibly  photography)  be  used  for  determining  y. 

Overall  Assessment  of  the  System 

Our  first  goal  with  the  remote-sensing-based  system 
was  consistency  of  estimates  among  areas  and  years.  We 
believe  that  this  goal  was  accomplished,  although  major 
differences  in  habitat  among  waterfowl  management  dis¬ 
tricts  demand  some  regional  modification  of  methods.  Our 
second  goal  was  maximum  use  of  existing  data.  We  be¬ 
lieve  that  this  objective  was  accomplished,  but  implemen¬ 
tation  of  the  remote-sensing-based  system  clearly  pointed 
out  that  certain  data  required  by  the  system  are  scarce  or 
lacking.  This  is  especially  true  for  reliable  estimates  of 
clutch  success  in  some  areas  and  in  some  habitats  (Shaffer 
and  Newton,  in  press).  Our  third  goal  was  rapid  execution 
of  complex  procedures  by  microcomputer  and  a  resulting 
database  that  documents  changes  in  the  estimated  parame¬ 
ters.  This  goal  was  partly  accomplished.  Because  the  re¬ 
mote-sensing-based  system  has  been  constantly  evolving, 
continued  modification  of  computer  programs  has  been 
necessary.  At  the  same  time,  technological  advances  in 
computer  development  and  remote-sensing  techniques 
have  moved  ahead  of  those  used  in  the  remote-sensing- 
based  system. 

Fluctuation — sometimes  violent — of  breeding-popula¬ 
tion  sizes  and  production  of  prairie  ducks  causes  problems 
for  waterfowl  managers.  The  manager  usually  attempts  to 
manage  such  things  with  the  amount  and  extent  of  harvest 
and  the  availability  and  quality  of  habitat.  The  success  or 
failure  of  management  is  usually  evaluated  in  terms  of 
duck  numbers  and  production.  Johnson  and  Shaffer  ( 1 987) 
demonstrated  the  difficulty  of  separating  weather  factors 
from  the  numerous  other  causes  of  fluctuation  in  duck 
population  sizes  and  the  need  for  long-term  surveys  with 
consistent  methodology.  The  remote-sensing-based  sys¬ 
tem  proposed  here  has  the  potential  to  provide  better 
evaluation  of  management  by  providing  data  that  are  es¬ 
sential  to  understanding  the  interaction  of  various  factors 
that  cause  fluctuation  in  duck  population  sizes. 
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As  the  Nation’s  principal  conservation 
agency,  the  Department  of  the  Interior  has  re¬ 
sponsibility  for  most  of  our  nationally  owned 
public  lands  and  natural  resources.  This  respon¬ 
sibility  includes  fostering  the  sound  use  of  our 
lands  and  water  resources;  protecting  our  fish, 
wildlife,  and  biological  diversity;  preserving  the 
environmental  and  cultural  values  of  our  na¬ 
tional  parks  and  historical  places;  and  providing 
for  the  enjoyment  of  life  through  outdoor  rec¬ 
reation.  The  Department  assesses  pur  energy 
and  mineral  resources  and  works  to  ensure  that 
their  development  is  in  the  best  interests  of  all 
our  people  by  encouraging  stewardship  and  citi¬ 
zen  participation  in  their  care.  The  Department 
also  has  a  major  responsibility  for  American 
Indian  reservation  communities.  ^ 


